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ABSTRACT 


The paper illustrates the use of continuum models in control design for stabilizing 
flexible structures. A 6-DOF anisotropic Timoshenko beam with discrete nodes where 
lumped masses or actuators are located provides a sufficiently rich model to be of inter- 
est for mathematical theory as well as practical application. We develop concepts and 
tools to help answer engineering questions without having to resort to ad hoc heuristic 
(“physical”) arguments or faith. In this sense the paper is more mathematically ori- 
ented than engineering papers and vice versa at the same time. For instance we make 
precise time-domain solutions using the theory of semigroups of operators rather than 
formal “inverse Laplace transforms ” We show that the modes arise as eigenvalues of 
the generator of the semigroup, which are then related to the eigenvalues of the stiff- 
ness operator. With the feedback control, the modes are no longer orthogonal and the 
question naturally arises as to whether there is still a modal expansion. Here we prove 
that the eigenfunctions yield a biorthogonal Riesz basis and indicate the corresponding 
expansion. We prove mathematically that the number of eigenvalues is nonfinite, based 
on the theory of zeros of entire functions. We make precise the notion of asymptotic 
modes and indicate how to calculate them. Although limited by space, we do consider 
the root locus problem and show for instance that the damping at first increases as 
the control gain increases but starts to decrease at a critical value, and goes to zero as 
the gain increases without bound. The undamped oscillatory modes remain oscillatory 
and the rigid-body modes go over into deadbeat modes. 

The Timoshenko model dynamics are translated into a canonical wave equation in 
a Hilbert space. The solution is shown to require the use of an “energy” norm which 
is no more than the total energy: potential plus kinetic. We show that, under an 
appropriate extension of the notion of controllability, rate feedback with a collocated 
sensor can stabilize the structure in the sense that all modes are damped and the energy 
decays to zero. An example, non-numeric, is worked out in some detail illustrating the 
concepts and theory developed. 
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1. INTRODUCTION 

The purpose of the paper is to illustrate the use of continuum models in 
control design for flexible structures: to provide the tools necessary to address relevant 
engineering issues. This is admittedly hazardous on two counts: on the one hand the 
complicated 3D geometry of realistic structures makes it almost impossible to use 
continuum models; while on the other hand the mathematics of continuum models is, 
for the most part, of mathematical interest only, and reduction to engineering practice 
is seldom undertaken. The alternative, universally the rule now, is to stay with the 
finite-dimensional “Finite Element” models. However the latter has the drawback 
that any control design is limited to specific numerical values of system parameters, 
and the dimensions, can be prohibitively large. 

The 6-DOF anisotropic 1-D Timoshenko beam model is a convenient compro- 
mise from both sides. As shown in [Noor and Anderson 1979, Noor and Russell 1986, 
Wang 1994, and Balakrishnan 1992] it is excellent for modeling lattice trusses. On 
the other hand the mathematical theory strikes a good balance between the trivial 
and the nontractable. 

The purpose of the control is to enhance the stability of the system, and the 
main interest in the theory centers on the modes and the damping attainable — 
the eigenvalue problem. A purely formal Laplace transform analysis can yield an 
entire function whose zeros in the complex plane are the eigenvalues. But the main 
difficulty is in determining the time-domain solution and the nature of the stability. 
Here is where it becomes necessary to use the theory of semigroups of operators and 
associated techniques from abstract (functional) analysis. In particular the system 
with control is no longer self-adjoint — the mode-shapes (eigenfunctions) are no longer 
orthogonal — and the problem of modal expansion of the solution must needs draw 
substantially on non-self-adjoint operator theory. A standard reference for the latter 
is [Gohberg-Krein 1969]. 

While some background in Hilbert space and liner r operator theory is assumed 
(and almost all that is needed is covered in [Achieser-Glassman 1966, Riesz-Nagy 1955, 
Balakrishnan 1981]) every effort has been made to make the treatment self-contained. 
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We begin in Section 2 by translating the basic beam equations into a canonical 
abstract “vibration” or “wave” equation in a Hilbert space. The main feature of 
the choice of state is the inclusion of the displacements at the discrete nodes (the 
“boundary” points) in it — as pioneered by the author for the Bernoulli beam model 
of SCOLE in [Balakrishnan 1991a]. Not unlike the FEM version, the “boundary” 
conditions for the elastic equations are so chosen (in particular making the stiffness 
operator self-adjoint and nonnegative definite) as to yield the correct form for the 

potential energy. 

Section 3 deals with spectral analysis. We show that the undamped structure 
modes are the zeros of an entire function and that the mode shapes yield an orthogonal 
basis for the space. We characterize the rigid-body modes showing that they span a 
six-dimensional space. We develop the Green’s function for the eigenvalue problem. 
We show the relation of the square root of the stiffness operator to the potential 
energy of the beam. 

Section 4 treats the time domain solutions of the structure dynamics, including 
“weak” and “strong” solutions. We introduce the notion of the “energy” inner product 
and show how the theory of semigroups of operators applies and the relation of the 

resolvent to the familiar Laplace transform. 

In Section 5 we show that, under an appropriate generalization of the con- 
trollability condition to infinite dimensions, rate feedback using a collocated sensor 
can stabilize the structure in the sense that all modes are damped and the elastic 
energy decays to zero. We show that the damping coefficient goes to zero as the 
mode number increases without bound. 

Section 6 is devoted to calculating the asymptotic modes of the structure with 
rate feedback. We construct the 6 x 6 matrix which defines the mode shapes, the 
modes being the zeros of the determinant. We show that the latter is an entire 
function of exponential type. We show that there are deadbeat modes, equal in 
number exactly to the (dimension) number of rigid-body modes. The oscillatory 
modes of the undamped structure remain oscillatory regardless of the control gain. 
Because of space limitation we do not go into detail on the root locus problem. Making 
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precise the notion of “asymptotic modes,” we show that the asymptotic modes are 
the “clamped” modes where all nodes are clamped — no displacement is allowed. 

Section 7 deals with modal expansion. Since with feedback control the modes 
sure no longer orthogonal we have to use a “biorthogonal” system. We show that the 
eigenfunctions do provide a Riesz basis and develop a “modal” expansion. 

Finally, in Section 8, we include a non-numerical example, albeit simple, to 
illustrate the concepts and theory developed in the paper. 


2. THE ANISOTROPIC TIMOSHENKO BEAM: 

DYNAMICS AND STATE SPACE FORMULATION 

We begin by describing the dynamics of a 6-DOF-1D Timoshenko beam artic- 
ulated with lumped masses, and/or control actuators at a finite number of “nodes” 
distributed along the beam, including the end points where the lumped masses may 
also be “offset.” This is a natural extension of the familiar 1DOF-1D Timoshenko 
beam models found in standard texts (e.g., [Timoshenko 1974], [Meirovich 1967]). 

The model does not include the inherent damping since we are concerned only 
with the damping attainable with control — we refer to this as the “undamped” 
structure. 

Next we choose an appropriate Hilbert space as the state space and define the 
mass-inertia operator, the stiffness operator and the control operator — extending 
the notions familiar in finite dimensions. In particular the definition of the stiffness 
operator is based on the (elastic) protential energy of the structure. With these 
definitions the partial differential equations translate into a “vibration” or “wave” 
equation in a Hilbert space which needs to be interpreted appropriately but provides 
the canonical model for the rest of the paper. 

The anisotroic Timoshenko beam model adopted here appears to have been 
introduced by [Noor and Anderson 1979], and refined in [Noor and Russell 1986], to 
model the lattice truss strctures used for deployment in space. Later [Wang 1994] 
showed how such models could be derived starting from the general elastic solid 
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equations using homogenization theory. Details of deriving the beam equations and 
the elastic and mass-inertia parameters therein illustrated by a specific example can 
be found in [Balakrishnan 1992], including a comparison with FEM for calculating the 
step response. To minimize complexity we only consider the case of a single flexible 
beam. An example of a multibeam structure modelled as interconnected Timoshenko 

beams is given in [Balakrishnan 1991b]. 

Let (x lt £ 2 , £ 3 ) denote the coordinates of a rectangular coordinate system, and 
let the beam axis be the £i-axis. We shall use “s” to denote the position along the 
beam: 0 < s < L, where L is the length of the beam. We consider the “uniform” 
case where the beam properties do not depend on s. Let Si denote discrete points 
(referred to as “nodes” ) along the axis 


Sj Sj+i, i 2, ■ • - , 771 1 

Si = 0; s m = L . 

Between nodes, that is to say for s t < s < Si+i, we have the basic anisotropic 
Timoshenko beam equations governing the beam displacements, u( - ), u(-), w( ). 

«(•) axial displacement (or, ii-component) 

u(-) bending displacement in the xi -£2 plane (or, X 2 -component) 
to(-) bending displacement in the £i -£3 plane (or, 13 -component) 

and the torsion angles <^>i(-) , ^3(‘) : 

(p^-) rotation angle about the aq-axis 
4 > 2 (-) rotation angle about the £ 2 -axis 
4>z(') rotation angle about the X 3 -axis 


are given by 


mult — Cult" — Ci 4 l>" — C\sw" — C1502 + Ci4<^3 = 0 
77122^ — C44U M — C14U" — €45^2 — C A5 W " = 0 
77133?!) — C55W" — C\zu" ~ C^ 4>'<2 ~ c 45 v ” — 0 
— C66<t>" — C36^>2 — c 26<^3 = ® 

77155^2 + m 564 l i + Cisti* + C 55 U/ — CZ6<t>\ + C5502 ~ c 33$2 — c 23<^3 — c 45<^3 = 0 
m66</>3 + "l56<^2 - Ci 4 u' + C44U' - C 2 6 <£i “ ^2302 + C44< ^ 3 _ ° 22 ^ ~ °^ 2 ~ 0 • 
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where the superdots denote time derivatives and the primes, the space derivatives 
(with respect to s). The matrices 


Cll 

C 14 

Cl5 


C66 

C36 

C26 

c 14 

C 44 

^45 

; C3 = 

c 36 

C33 

c 23 

Cl5 

C 45 

c 55 


C26 

C23 

C 22 


are both strictly positive definite. We shall also use the notation: 



0 — C15 C14 


0 0 0 

c 2 = 

0 — C45 C44 

0 - c 55 C45 

; C 4 = 

0 C55 — c 45 

0 043 C44 


With / denoting the 6x1 (column) vector: 

u 

v 

f= W 

1 4>i ' 

4> 2 

4> 3 

these equations can be conveniently rewritten in the vector form: 


Mof — A-it" + Aif' + Aof = 0, Si < s < Sj + i 


where Mq is the mass/inertia matrix 


M o 


mu 

0 

0 

0 

0 

0 


0 

77122 

0 

0 

0 

0 


0 

0 

77I33 

0 

0 

0 


0 

0 

0 

m 44 

0 

0 


0 

0 

0 

0 

m S 5 

m 5 s 


0 

0 

0 

0 

mss 

mss 


where M is also required to be strictly positive definite; 


and where / stands for 


Ai 


Ci 0 

0 Cz 


Ai 


0 C 2 

— C2* 0 


^0 


0 0 

0 c 4 


f(t, s), 0 < t, 0 < s < L. 


(2.2) 


( 2 . 3 ) 
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Here the m tJ are the “mass inertia” coefficients and { c x] } the “flexibility coefficients. 
They are constrained by the positive definite requirements of C\ and C3; in particular 
ca > 0, mu > 0 for all i. 

The nodes {sj are points where the slope — spatial derivative /'(•) — is 
discontinuous because of controllers or lumped masses located there. The displace 
ments at the nodes thus have to be included as part of the definition of the state. 
Let I 2 (0, L) 6 denote the I 2 -space of 6 x 1 vector functions /(•) and let H denote the 
Hilbert Space 

H = I 2 (0, L ) 6 x (R 6 ) m 

where we use the notation (f? 6 ) m rather than fl 6 ”* to indicate 6 x 1 vectors replicated 
m times. For elements x in Ti. we shall use the notation 

Hi| 

where / € L 2 (0, L) 6 and b 6 (R e ) m - To avoid possible confusion let us note explicitly 
that the inner product in H is given by 

m 

[i, y] = [/, 0] L2(o, l ) 6 + 


where 




y = 


9 

c 


C 1 
C 2 

Cm 


(To avoid excessive notation we may often delete the signature under the inner prod- 
ucts in the sequel if they are clear from the context.) The elements of H will be our 
“states.” We begin with the stiffness operator since the mass operator will depend 
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on the control masses, the stiffness operator defined as the “differential” operator A 
with domains, denoted T>(A), in H given by: 


V(A) = 


x = 


where defining /j(s) = /(s), s* < s < s i+l) i = 1 , . . . ,m — 1 

fi(-), f^) 8md f l' (-)eL 2 ( Si ,s i+1 ) 6 


lim fi(s ) = lim f i+1 (s ), * = 1, . . . , m - 1 


b = 


f(s l) 


/(Sm) 


In other words the functions /(•) in V{A) are “piecewise smooth”; continuous in 
0 < s < L, but the derivative can have a jump discontinuity at each node S{. We 
refer to 6 as the “boundary value.” The operator A is defined by: 


Ax = y\ x = 


f 

b 


y = 


(2.5) 


where 


g(s) = -A 2 f”{s) + Aif'{s) + A 0 f(s), s t < s < s i+i , i = 1 , . . . ,m-l. (2.6) 


and c is defined by: 


c = 


~LJ( 0+) - A 2 f'(0+) 
-M(f(s2+) - /'($ 2 -)) 

~ Mif'iSm- 1 +) ~ )) 

Lif{L-) + A 2 f'{L-) 


(2.7) 


where 


L i — 


0 -C 2 

0 0 


We shall find it convenient to use the notation: 


c = A b f. 


This definition is made so that we get the right expression for the potential energy 
— so that 


[. Ax , x 


= Potential Energy 


( 2 . 8 ) 
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- 

u' 


u' 

“ 

rL 


<*>; 




Cl 

u' - <£3 

> 

v' ~ 03 


ds ■ f / 

1 A 

c 3 

02 

> 

02 

<^3 



w' + 02 


It?' + 02 


JO 


<^3 


- 


ds. 


(2.9) 


where the potential energy of the structure is defined by 

s: 

Let us verify (2.8). Integration by parts yields: 
fS ' + W"( 5 ),/( S )] ds 

= [^/'(St+l-^/^t+l-)] - [^ 2 /'( 5 »+)./( s *+)] _ j s [Mf'(s),f (s)] ds. 


Hence 

rL 


- f \A 2 f"(s)J(s)]ds 
Jo 

m— 1 m— 1 pL 

= -S>»/Wi-)./(*+i-)l + I> 2 /'(*+). /(*+)] + J 0 /(»)]<**• 


But 


And using 


m — 1 

- EM2/W1-). /(».«-)] = 

1= l 


*=2 


f(Si-) = f(Si+ ) = f(Si) 


( x being in T>(A)), we have that 

m— 1 m-l 

~ S [A2/'(s*+1-)i f(Si+l~)] + [^2/'(Si+). /(«i+)] 

i=0 t=l 


m— 1 

= [A 2 f'(L-),f(L-)} + j2[i4 2 (/W)-/'(*i-))./(*i)][^2/'(0+),/(0+)]. 

2 

Also, noticing that 

Ai — L[ — Li, 


we can calculate that 


Vi/'M./MM 8 = -[i 1 /(i-)./(i-)] + (ii/(0+),/(0+)] 

+ (V/M, /'«]<!»+ Am/'m./mw*- 

JO JO 
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But 


[Abf,b] = l(Lif(L-) + A 2 f'(L-))J(L-)} 

m - 1 

+ 51 [^(/'(s.-) - /'($»+)), /(«»)] 
2 

-[Ii/(0+) + d 2 /'(0+),/(0+)]. 


Hence 


[■Ax.i] — [Abf , 6] + [(Aof + A\f' — Atf"), f(s)] ds 

= jfW'(s), /'(«)] * + jfW(«), /'(*)] 

+ ^ l L if'( s ), f(s)] ds + f Q [^o/(s),/(s)] ds. 


Hence 


where 



H 


Ci 0 0 -Ci 

0 C 3 0 0 

0 0 M 

—Cz* 0 


from which (2.9) readily follows, noting that 


( 2 . 10 ) 


C2 = —C\Dz 


D 3 


0 0 0 

0 0-1 

0 1 0 


C 4 = DlCiDz 


or, equivalently: 


A 0 = L\AZ l h. 


This technique of defining the “boundary conditions” on the “differential” operator 
appropriately to make it self-adjoint and nonnegative definite should be compared 
with the FEM method where (in the Hamiltonian) the potential energy is specified 
and the stiffness matrix derived therefrom. 
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Remark 

By our definition, a node is a point where there is a controller or lumped mass. 
There may be neither at the end points s = 0, s — L, and if so, they are not included 
in the nodes, and additional conditions need to be imposed. Typically, either 

/( 0) = 0 (clamped at zero) (2.11) 

or 

L x /(0) + A 2 f'(0 ) = 0 (free at zero) (2.12) 

and similarly at s = L. 


Mass/Inertia Operator 

To obtain the mass/inertia operator we have to specify the control mass/ 
inertia and end masses, possibly offset (modeling antennas). Let us begin with the 
interior controllers, force (reaction jets, proof-mass actuators) and moment (cmg’s). 
Then we have for i / 1 or m 


M b J{t,s i ) + A 2 {f'{t,s i -)-f , {t,s i +)) + U i (t) = 0 (2.13) 


where 


M m = 


m* 0 0 0 0 0 

0 mi 0 0 0 0 

0 0 mi 0 0 0 

0 0 0 

0 0 0 Ii 

0 0 0 


and Ui is the 6x1 vector of force and moment controls at s = s t ; m, is the control 
mass and /* the corresponding moment of inertia (matrix). For the ends, allowing for 


offset end masses mo at s = 0 and mi at s = si and Iq and I c the moments of inertia 


of antennas and controller, each about the center of gravity, and r 0 and r L the 3x1 
position vectors of the centers of gravity respectively of the end masses we have at 


s = 0: 


0+) — Lif (t, 0+) — A 2 f'(t, 0+) + Uo(t) = 0 


(2.14) 
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where 


and 



and where m\ is the controller mass; and finally Uq is the control vector (forces and 


moments) . 


Similarly for the end s = L, we have 


where 


M blL f(t,L-) + L 1 f(t,L-) + A 2 f'(t,L-) + U L (t) = 0 

L 

r L ~ r L ,2 

r L,Z 


(2.15) 
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and m m controller mass and Ui is the control vector of forces and moments. Again, 
(mi + m m ) is the total moving mass. 

Let Mb denote the composite matrix of all the control and end masses/ inertia: 


M b 


Mb, o 

Mb, 2 


M b,m— 1 

Mb,L 


The mass-inertia operator, denoted M, is then defined by 


M x = y; x 


f 

b 


y = 


M 0 f 

M b b 


( 2 . 16 ) 


and is a linear bounded self-adjoint nonnegative definite opeator H onto H with a 
linear bound inverse: 

M~ l x = 


and we note that 


Also 


Mo l b 

M; l b 


M b l b = 


Mbflh 


\/~M x = 


M b y*>i 

Mb,lbm 

y/M~of 

y/W h b 


\J~Mb b = 


\J M b ,i b{ 
\J M b,L b m 


Control Operator 

The control operator maps the control inputs into H. We denote it by B. We 
note that it is possible that not all nodes may have controllers. Let the number of 
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controllers — force, moment or proof-mass — be denoted m c , 


m c < 6 m . 

Then we can regard any control u as an m c x 1 column vector. We define 

B u u = b 


where 


b = 


B\u 

B m u 


where each Bi is a 6 x m c matrix such that B*Bi is nonsingular and further 



* ^ 3 


or 


BiB*= 0, i?j. 


In particular it follows 




We shall use the notation 


BiBfa 

BmB m b m 


BiB * = Di 


where the D j, which are 6x6 each, may be taken to be diagonal with entries 1 or 0. 
Finally 


Note that if 


we must have: 


BB* 


f 


0 

b 


B«B' u b 


B*x = 0, 


x = 


/ 

b 


B* u b = 0 
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or b must be such that: 


B-bi = 0, i = 1 , . . . ,m, 

or 

DA = 0, t = 1, . . . , m. 


Another equally useful representation for BB* can be obtained in the following 
way. Let {ej, t = 1, . . . ,m c , denote the unit coordinate vectors in TZ mc . Then let 

B(i) = Bei . 


Then 


[S(l). B U )\ = 0, ■ f J 


BB- = 


where the B (*) are orthogonal, and each B^) is nonzero. 

With these definitions we can assemble the canonical “state space” version of 
the dynamic equations (2.1), (2.13), (2.14), (2.15): 


Mx(t) + Ax(t) + Bu(t) = 0, t > 0 


(2.17) 


where, formally, we have taken: 


x{t) = 


/(«, •) 

b(t,-) 


We recognize (2.17) as “an abstract wave equation in a Hilbert space”; the precise 
relationship between the original space-time dynamic equations and this abstract 


version will be clarified later, in stages. 


3. SPECTRAL ANALYSIS: MODES/MODAL EXPANSION 

Before we proceed to “solve” (2.17) we need to examine the spectrum of the 
stiffness operator A. The eigenvalues of A (with respect to the mass operator M) 
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are the modes of the undamped structure. A crucial result is that the modes are the 
zeros of an entire function — the determinant of the m x m “condensed dynamic 
stiffness matrix.” The corresponding “mode shapes,” the eigenfunctions, yield an 
M-orthogonal basis for the Hilbert space. We also characterize the “rigid body” 
modes, the eigenfunctions corresponding to the zero eigenvalue. We also obtain the 
“Green’s function” for the eigenvalue problem. Finally we define the square root of 
the stiffness operator and indicate its relation to the potential energy. 

Elementary Properties of the Stiffness Operator 

Let us begin cataloging the elementary yet crucial properties of the stiffness 
operator A. We have already seen that it is self-adjoint: 

[Ax, y] = [x, Ay] for x, y in V(A) 

and nonnegative definite 


[Ax, x] > 0, for x e T>(A). 

The next property is that the domain of A is dense in Ti. In other words, while not 
every element x in H is in the domain of A as the definition of A clearly shows, we 
can find elements in the domain of A that approximate x as closely as needed. More 
precisely, we can find a sequence {x n } in T>(A) such that 

||x x n || * 0. 

This is a feature of “differential” operators; but since in our case we also have “bound- 
ary values” to contend with, we shall present a formal argument. Thus let 

/ . 

6 

Let 

y*(^) /(^)t s* ^ s ^ 5j+i- 

It should be noted that b is arbitrary and is not necessarily the vector of “boundary” 
values of /(•) at s = s»; the latter need not of course be even defined. But for 
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each i, we can find a sequence of functions {/t,n( - )} suc h that /t,n( s ) is infinitely 
differentiable in Sj < s < Sj+i, and 


and 


r +l Wfiis) ~ /i,n(s)|| 2 ds — 0 as n — ► oo 
J Si 


fi,n{ s i ) — /i,n( s t+l) — ^i+1- 


Defining 


/ti(s) — /t,n(®)i < S < 5t+l, 

fn(Si) ~ /i,n(®)i 


we see that /„(•) is continuous and 



belongs to the domain of A, and 

| |x — arn|| 2 — 0 as n — ► oo 


as required. 


Eigenvalues and Eigenfunctions 

If we proceed formally and take the Laplace transform of (2.17), setting 

<p = [ e~ xt x(t) dt , Re A > 0, 

Jo 


we obtain for fixed A: 

A 2 M <t> + A4> — ip = 0 


(3.1) 


where ip is an element of H . The solution of this equation plays an important role in 
the theory. We now formulate this more precisely as: find <p in V(A) such that 


X 2 M<p + A<p = ip, ip € H 


(3.2) 


for given A. We see that in order for (3.2) to have a unique solution the “homogeneous” 
equation 

A ^M<p-\-A<p = 0 (3-3) 
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can only have the zero solution. Rewriting this equation as 


A<p = -X 2 M<p 


we may consider more generally 


A<t> = 7 M<f>, 4> / 0 . ( 3 . 4 ) 

Here 7 is called an eigenvalue of A with respect to the mass matrix M, and 4> a 
corresponding eigenfunction. FYom 


0 <[A</>, <f>] =i[M<t>,<t>) 


it follows that 


7 > 0. 


If 7i) 7o are two distinct eigenvalues with corresponding eigenvectors <pi, fa respec- 
tively, we have 


\A<t> 1, <t> 2] = fa] = [fa, A<t>^ = 7 2 [M^>i, <f> 2 ] 

and hence we must have: 

[■ M<f>i , <fo] = 0 

or, eigenfunctions corresponding to distinct eigenvalues must be M-orthogonal. 

We shall show now that the set of eigenvalues is nonfinite, countable and can 
be taken as 

{wit}, < Uk+i, k = 0,1,2,... 

and 

Uk — » 00 as k — ► 00 


and 


A<t>k = u/fcM 4>k- 


First let us consider the eigenvalue zero. 
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Rigid Body Modes 

Because we are considering the case where both s = 0 and s = L are nodes, 
we shall show that zero is an eigenvalue. An eigenvector corresponding to the zero 
eigenvalue is called a rigid body mode. This is because if 

A4> = 0 

we have 

[A<f>, </>} = 0 

or, from (2.9) the associated potential energy is zero. Let 



where b is the boundary vector, 4> being in V(A)\ if the potential energy is zero, we 
have from (2.9) that /(•) must be of the form: 

u(s) = u(0) 
v(s) = v(0) + 5^3(0) 
tu(s) = U>(0) — S02(O) 

*i(«) = *i(0) 

0 2 (s) = <fo(0) 

4>3(s) = </> 3 ( 0 ) 

The derivative f'(s) is continuous in 0 < s < L. Note that the “free-free” boundary 
conditions (2.12) are satisfied at both ends, and /(•) is the same regardless of the 
number and location of the interior nodes. Also the dimension of the eigenfunction 
space is six. In our notation we shall set 

u>o = 0. 

Nonzero Eigenvalues 

Let us now consider the general case in the form 

A 2 M<t> + A<p = 0. 
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Let 



Then 


A 2 M 0 / +5 = 0 


where g(-) is defined by (2.6), and 


A 2 M b b + A b f = 0. 


To solve (3.3), let 


A(X) 


Then 


0 h 

AI\Aq + \ 2 Mq) A^Ax 


m 




/(Si) 
f'( s i +) 


Si ^ S ^ ^i+1) 


1,2,. . . ,m — l. 


Let 


where 


a 


f(s l) 

ms i) 

m* 2 ) 


Sf'(Sm-l) 


m*i ) = /'(*+) - /w-) 


and we set 

. /'( o-) = 0. 


Then 


f(s) 


/(0) 

+ V e^ (A)(s " Sj) 

0 

f'(s) 


0 


«/'(*,) 


Si < s < s i+ 1 , i = l,... ,m — 1. 


Hence we can write 

6 = L(A) a 


(3.6) 


(3.7) 


(3.8) 


(3.9) 


(3.10) 
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where L{ A) is a 6m x 6m block lower-triangular matrix: 

L{ A) = {4j(A)} 
tij( A) = 0, j > i 
tnM = h- 

For i > 1 

4,1 (A) = Pn(A,Si) 
4.2(A) = Pi2(A, Si) 
4.3 (A) = Pi 2 (A, Si - s 2 ) 


where 


Note that 


4,t-i(A) = Pi2(A, Si - Si-2) 
4 ,»( A) = Pi 2 (A, Si — Sj-i) 




Pn(A, s) Pi 2 (A,5) 
P 2 l(A,s) P22(A, s) 


3 3 

P 2 i(A,s) = — Pu(A,s), P 2 2(A,s) = ^Pi2(A, S ). 


det 


L(A) = detPi 2 (A,s 2 ) ••• det P l2 (A, s* - s^) ••• detP l2 (A, L — s m -i)- (3.11) 


Also, we can express A b f in terms of &, using (3.10). Thus we can write 


A b f = K(x)a 


where K( A) is also a 6m x 6m (block-lower-triangular) matrix and hence (3.3) yields: 

(a 2 M 6 L(A) + K(\j) a = 0. (3.12) 

Hence the eigenvalues are the roots of 

det(A 2 M fc L(A) + K( A)) = 0. (3.13) 

The left side defines an entire function of the complex number A and can have at 
most a countable number of zeros which must grow without bound in magnitude, if 
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nonfinite in number. The corresponding eigenfunction is determined by (3.10). Also 
(3.12) shows that the dimension of each eigenfunction space is finite, not more than 

6 777. 

We can now deduce that the set of eigenvalues which is countable is not finite. 
For if it were finite, the number of eigenfunctions would be finite. But the eigenfunc- 
tions must be complete in the M-inner product. If not, there must be a function ip 
which is Af-orthogonal to all the eigenfunctions. The class of such functions form a 
Hilbert space and A must map a dense domain of this space into itself. But A being 
self-adjoint, we can now recall the fundamental result (see, for example [Riesz-Nagy 
1955]) that is must have at least one eigenvalue, which leads to a contradiction. But 
the dimension of H 2 being nonfinite, we see that the set of eigenvalues cannot be 
finite. In particular we see that the eigenfunctions are complete, and M-orthogonal. 
We can therefore make them an M-orthonormal basis for H. For each eigenvalue u%, 
let Pk denote the projection (operator) corresponding to the eigenfunction space: 

A(P k * ) = v&Pk*). 


Since the dimension of each eigenfunction space is finite it is convenient to continue to 
use {^fc} to denote the M-orthonormalized basis in each eigenfunction space, “count- 
ing each eigenfunction as many times as the dimension,” as is customary. Then we 
. can write 

OO 

x = (3.14) 

0 

where 

Q k = [x,M<p k ] 

and 

£bt| 2 = [Mx,x]. 

0 

Hence 

OO 

M~ l x = ^2[x,<pk]<Pk 
0 

or, 

OO 

1 = 'Y![x,<PkW <p k - (3.14a) 

0 
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Using (3.14) and (3.14a) we have 

[x,x] = ^2\x,4>k][M(f>k,x] 
0 


conversely 

OO 


X = ^Ctic&k € H 
l 

if and only if 

OO 

£l<**| 2 < 00 
1 

and 



[x, Mx] = X^l afc | 2 - 
1 

Since 

A4>k = u}\M<f>k 


we have that 


it follows that 



x € V(A) 


if and only if 

OO 

£|[x,M^]| 2 ^ < oo. 

l 

Also for x, y in H: 

OO 

[x, My] = ^[x, y] (3.15) 

o 

[x,y] = x,<t>k\[M<i>k,y]- (3.15a) 

o 


What we are exploiting here is the fact that the sequences {<£*}, {M <p k } are “biorthog- 
onal” and complete. We shall return to this concept later in Section 6. We refer to 
(3.14) as a “modal” expansion. Expanding on (3.14), let 


23 



Then (3.14) yields: 


Taking 


we have 


Hence 


and hence it follows that 


and hence that 


In particular therefore 


b 


OO 


^2\x,M<pk]bk' 

o 



0 

b 


b = ^[b, M b b k ]b k . 
0 


O 6 * M ^)i 2 < 00 


OO 




< 


OO 


00 


ElMI 2 < 


OO . 


bk —* 0 as k — * 00 . 


Returning now to (3.10) we have 

i 

f(s ) = Pn(A,s)/(0) + ]r p i 2 (A, s _ 8j )6f'(8j), Si<s<s i+ 1 

j=X 

which we can express as a linear transformation: 


/ = £(A)a, 


where £(A) maps ( R 6 ) m into £ 2 ( 0 , £) 6 - Let 


Then 


h = £( A)a lf a l€ (R 6 r 
/ 2 = £(A)02, 02 G (R 6 ) n - 


i>\ = 


£(A)Oi 
L(A) Cbi 


i> 2 


£(a)02 

L[\)(h 
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(3.16) 


(3.17) 



are elements in T>(A) and 


[Ail) 1,^2] = [Abfi,h] - A 2 [M 0 /i, / 2] 

= [KiX)^, L{\)(h} - A 2 [M 0 /i, h] 

= tyuAfa] 

= [L{\)a u k(x)02] - a 2 [m 0 /i, hi 

Hence for A such that A 2 is real 

L(\yK(\) 

is self-adjoint and is nonnegative definite if A is real, since 

0 < [Aipi,Tpi\. 

Also L( A) is an entire function of the complex variable A, and 

det|L(A)| = 0 


(equivalently 

detP 12 (A; Sj - Sj_i) = 0 

for some i > 2, from (3.11) ) for at most a countable set of A. And, omitting this 
set (which are recognized as “clamped” modes of the structure, where every node is 
clamped so that the displacement is zero), we have 

a = L{ \)~ l b 


and 

(a 2 M 6 L(A) + A(A)) a = (x^Mb + K (A)L(A) -1 ) b. 
Hence we can express (3.12) as: 

(A 2 M 6 + A(A)L(A) _1 )6 = 0. 

Following [Wittrick- Williams 1971], we shall call 

A 2 M b + K (A)L(A) -1 
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the “condensed dynamic stiffness matrix.” Let 


T(X) = K(X)L(X)- 1 . 


Then we have: 

[A4\M = [AbfiM ~ A 2 [M 0 /i, / 2 ] 

= [K( A)L(A)- 1 i. ll b 2 ] - A 2 [Mo/i, h) 

= [4> i, A<t> 2 ) 

= [d lf /f(A)L(A)- 1 6 2 ] - A 2 [M 0 /i ( f 2 }. 

It follows that for A such that A 2 is real, T(A) is self-adjoint; and if A is real, also 
nonnegative definite. Also, if 


K(x)a=o, a^o, 


taking 


we have 


f = c{ A) a, b = L{x)a, 



A b f = o 


and hence 


\A<t>, 4>\ = [A b f,b] - X 2 [M 0 fJ] 

= -A 2 [M 0 /,/] 

which is impossible unless A is pure imaginary. Zero is an eigenvalue of K (0). By the 
“zeros of the dynamic stiffness matrix” we mean the roots of 

det(A 2 M 6 + T(A)) = 0 

which are of course the same as those of 

det (A 2 M 6 L(A) + A (A)) = 0. 
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The poles of the condensed dynamic stiffness matrix are the zeros of 

det L(A), 

or the “clamped” modes of the structure. The dimension of the eigenfunction space for 
A ^ 0 is equal to 1 if the condensed dynamic stiffness matrix has distinct eigenvalues. 


Green’s Function 

We are now ready to solve the non- homogeneous equation: 

(A 2 M + A)4> = it. 


(3.18) 


In fact for A such that 


det (a 2 M 6 + T(A)) ^ 0 


or 


A 2 * -u& 

we can calculate the solution in two ways. First we can obtain a “modal” solution 
using the modal expansion (3.14). For this purpose we note first that: 

<t>k 


(\ 2 M + A)- l M<f> k = 


A2 + 


OJt 


since 


Hence we obtain: 


(A 2 M + A)4>k = (a 2 + ul) M<t>k. 


and 


(A 2 M + A)~ l x = 53 


OO 

[(A 2 M + A1) _1 X, y] = 53 


[x, 4>k\<f>k 
A 2 + 


OJt 


[x,4>k][4 >k, y] 
A 2 + 


ut 


2 


For A 2 positive we see that 

(A 2 M + A)' 1 

is self-adjoint and nonnegative definite, 


[(A 2 M + A)- 1 4k, M<Pk} 


1 


A 2 -I- 
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(3.19) 


so that 

oo oc -r 

+ M* k \ = z^T w . 

We shall show that the right side is finite. For this purpose we need to derive the 
Green’s function for the eigenvalue problem (3.18). 

For this purpose let 


<P = 


ip = 


Then (3.18) becomes 


A 2 M 0 / + h = g 


where h(-) is defined by (2.6), and 

X 2 M b b + A b f = c. 

With v4.(A), as before we see that we can solve (2.38) as 


(3.20) 


f(s) 

= g A(\)S 

/( 0) 


,° 1 + 

P(^) 

/'(*) 


0 

j= 1 

If'isj) \ Jo 

0 

( 


da. 


Or, with CL, C( A) as before, we have: 

/ = C{ A) a + Af(X )g 

where the operator M(X) mapping I/2[0, L] 6 into itself is defined by: 

h=Af{\)g\ h(s) = [ P n (X,s — a) g(a) da, 0 < s < L. 

Jo 

Also, we have 


where 


b = L{X)a + N(X)g 

MO) 

HL) 


N(X)g = 


and 


A b f =K{X)a+M{X)g 


28 



where 


M{X)g = A b (Af{X)g). 


Hence (3.20) becomes 

(A 2 M b L{X) ■+- K( A)) a = c - N (X)g - M (X)g 


and hence for 

A 2 + 7^ 0, 


we have: 

a = (A 2 M fr L(A) + if(A))' 1 (c- AT(A) ff -M(A)p). 

Hence finally 

/ = £(A)a(A) + ^(A)*/ (3.22) 

6 = L(X)a(X) = (A 2 M fe + r(A))) 1 (c-N(A)p-M(A)ff) (3.23) 

where 

a(A) = (A 2 M b L(A) + KiX))- 1 (c- N{X)g- M{X)g). 

From (3.22) we see that we can write 

r L 

f(s) = / G(X, s, a) g(a) da, 0 <s<L 
Jo 

b = (a 2 M 6 -I-T(A))” 1 c 4- G b {X,a)g(a) da = N(X)f 

where the kernel G(X,s ) a ) is continuous in 0 < s, a < L and G b (X, a) is continuous 
in 0 < a < L. It follows from this 

(A 2 M + A)' 1 


is nuclear (see [Balakrishnan 1981] for the definition) and in particular (3.19) is finite 
and hence it follows also that 


OO 1 

< »■ 

1 u k 

Of particular interest to us is the solution of (3.18) when rp is of the form 


ip = Bu 
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and (3.22), (3.23) simplify to: 


/ = C(X) (A 2 AfyL(A) + tf(A)) 1 B u u (3.24) 

b = (A 2 M 6 + r(A))" 1 B u «. (3.25) 

Hence 

B'{\ 2 M + A)~ 1 B = B;(A 2 M & + r(A)) -1 B„. (3.26) 


The Square Root of the Stiffness Operator 

Finally we note that we can define VI as a self-adjoint nonnegative operator. 
This is treated in standard texts — e.g., [Riesz-Nagy 1955]. It is known that 

D (VI ) D D(A). 

Unfortunately we cannot use the expansion (3.14), since we cannot evaluate 

VI 4>k or Va M<t>k- 


However we do know that 


V~A<fik, VI <t>j 


[A<p k , <f>j ] 

0, k ^ j 
Vfc, k = 3- 


Using this we have 


which if x € 'D(A), is 


OO 

/4i, VI x = y^u 2 |[x, 


(3.27) 


= [Ax, x\ — 2 (Potential Ener.jy). 


Hence we see that if 

x ev (VI) 
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then we can define 


(Potential Energy) 


llv'I j|1 2 

2 


Thus we can extend the definition of Potential Energy for all x in T> (\/A ) , which is 
larger than T>(A). 

Since the stiffness operator A has nothing to do with the mass operator M , 
we may consider the case 


M = Identity 


in which case, our eigenfunctions become 

A<pk = & k 4> k- 


Then in terms of these eigenfunctions we can define 

o° 2 

V (v^) = z £^|[s,^]| < oo (3.28) 

and 

00 r •, 

■/a x = ^ u>k |z, 4>k\ <f>k 

1 

[Ax,x] = f^|[z,^*]| 2 

which is then also the domain on which the Potential Energy can be defined. However 


we need a characterization similar to that for A , which in particular does not invoke 
eigenfunction. For more on this see [Balaknshnan 1990]. 


4. TIME-DOMAIN ANALYSIS 

We have tools enough to consider the “time domain” solution of the dynamic 
equations: (2.14) and in turn (2.1). Unlike the finite-dimensional case, we can have 
more than one kind of solution depending on the interpretation of the equation (2.14) 
and the properties demanded of the solution. For the most satisfactory form of 
solution, we need to introduce the notion of the “energy norm” space in which the 
norm is determined in terms of the total system energy — kinetic-plus-potential — 
and the theory of semigroups of operators applies. 
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We begin with the weaker notion first. 


Weak Solution/Modal Solution 

With respect to (2.14) we have to specify first how the time-derivative therein 
is to be defined and whether the solution x(t), 0 < t, which will depend only on the 
initial conditions at t = 0 and the input u(-), is required to be in the domain of A. 
Whatever the definition, we must have that for every <p in T>(A): 

[Mx(t), <}>] -I- [Ax(t), <fr\ + [Bu(t), 4>\ = 0. 

We can rewrite this as 

d 2 

— [x{t),M4>] -i- [x(i), A<t>) -I- [Bu(t), <t > ] = 0. (4.1) 

As for the input u(-), we assume that 

f |u(t)| 2 dt < oo, for every T, 0 < T < oo . 

Jo 

That is, u(-) 6 L 2 [0, T] mc for every T < oo. Note that in (4.1) we only need the “weak 
derivative” (see [Balakrishnan 1981] for definition) and further we have circumvented 
the requirement that x(t) € TJ>(A). By a “weak solution” of (2.14) we mean a function 
x(t ) , t > 0 such that for every <p i n ^(A): 

[x(t) t M<t>] 

and its derivative are absolutely continuous in t > 0, satisfies (4.1) and the initial 
conditions: 

[x(t), M <p] —>■ [xi, M d>] ast-+ 0+ 

-7- [x(t), M<t>] — > [x 2) M<(>\ as t — * 04- 

at 

where xj, x 2 are given elements in H. 

Even weaker (and easiest to construct) is the notion of a “modal solution” 
where we only require that (4.1) hold for every eigenfv nction 4>k ■ In this case (4.1) 

becomes 

d 2 



Let 

o/t(t) = [x{t),M<t>k\ 
u fc (t) = [u( 

Then for each fc, we have the ordinary differential equation 

ajt(f) + w 2 a fc (t) + u k {t) = 0, a.e. t > 0 (4.3) 


with the initial conditions: 

flfc(O) = [xi ,M<t> k ] 

dfc(0) = [i2, M<f> k ] 

which we can solve, to yield, for u k ^ 0: 

ak (t) = Q fc( Q ) s ^ na ^ + a k (0)cosuj k t + jf d- sin u k (t - s) u k (s) ds (4.4) 
and for u k = 0: 

a k {t) = a k { 0) + MO) + J\t - s)u k (s) ds. (4.5) 


It is easy to verify that 


OO 

£|a fc (t)| 2 < oo, 0 < t < oo . 

0 

Hence we can define 

x(t ) = 5Zafc(t)<^ifc. (4.6) 

o 

Thus defined, we do have 

Ofc(t) = [x(t), M <p k ] 

and *(•) satisfies (4.2). Hence we have a “modal” solution. However *(•) need not 
qualify as a weak solution in general and additional restrictions will need to be placed 
on the initial conditions — on Xi,i 2 - Thus for 


u(t) = 0, t > 0 


we would need to require that 


x\ e V{A) 
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so that 


In this case 


and hence for 4> in V{A) 


^a fc (0) 2 <4 < oo 


£M*)I 2 + £m*)I 2 < °° 


•^[x(t),M<t>] = Y,a k {t)[<t>k, M <fi] 


°° 1 
= 5Z ajk ^) (- u k)M<f>k, 4>\ 


= -^,ak(t)[4>k, A<f)\ 

0 

= -[x(f), A4>] 


as required. 

A better technique, avoiding these special considerations, is to go over to “state 
space” solution. 


State Space Solution: Need for Energy Norm 

To proceed further with (2.14), let us cast it in “state-space” form. Thus we 


let 


Y(t) = 


x(t) 

x(t) 


and, formally, (2.14) goes over into 


Y(t) = AY(t)+Bu(t) 


(4.7) 


where 

Y(t) eHxH, 


0 7 

-M~ l A 0 

0 

-M^B 
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t/i yi € 'D(A) 

V(A) = Y = V , V ~j ■ 

t/2 J/2 ^ /l 

We interpret (4.7) in the weak sense, satisfying also the initial condition: 


||Y(t) - Y\\ — 0 as t — ♦ 0 + . 

In addition we require that for each t, Y (t) is continuous with respect to Y : 

11^(011 — 0 as||Y||-0. 

It is shown in [Balakrishnan-Triggiani 1993] that this is impossible unless we change 
the space Ti x Ti — change the product-norm, which is physically meaningless, to the 
“energy” norm. The total energy 

= Potential (Elastic)Energy + Kinetic Energy 

[Ay i, yi] [My 2 , y2] 

”22 


where 


Y = , yieV(A). 

y 2 


As we have seen, we can extend the definition of Potential Energy to *D (y\/~A ^ . Also 


the total energy vanishes for 


Y = 


, yi € Nullspace of A. 


Let Tii denote the M-orthogonal complement of the null space of A: 


Tii = [x I [Mx,<t>] = 0, 4> e Nullspace of A}. 


Then 


<t>k ^ f° r every k , u k A 0. 


In other words we consider Ti under the M - inner product, 


[x,v]m = [Mx,y] 
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which is equivalent to the original inner product since M has a bounded inverse. We 
now define the energy space He by 

H e = (Z>(\ /A)nHi) x H 

with inner product defined by 

[Y, Z\e = [VAyi, '/A zij + [Afy 2l 22] 


where 


Note that 


Y = 

Vi 



; z — 


V2 


z 2 


[Y, Y]e = 2 [Total Energy] 

and hence the name “energy norm” space. We shall show now that He is actually 
complete. Let 



be a Cauchy sequence in He- Then 


||Vn - Y m \\% = |v7(x n - x m )|| 2 + ||\/m (z n - z m )f . 
Since \/~M has a bounded inverse, 


z n — > z in H. 


Now 

Hi, 

and \Ta restricted to Hi has a bounded inverse. Hence 

x n = (\/I) 1 '/A x n 

converges, and denoting the limit by x, we have 

x = (\/A ) (hm\/]4x n ) 
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and 


xeV (VI) DHv 


Hence He is complete. 

Let ^ — 1, - • • , 6, denote an M-orthonormal basis for the null space of A 

[M(po,k, 4>Qj\ — Sj- 

Define the projection operator (self-adjoint in the M-inner product) by 

6 

Vox = ^[x, M<t>o,k]<t>o,k 

l 

t 

and let 

V 1 = I - Vo 


so that 

[Vix, M4 > ] = 0, if V 0 4> = <f>- 


Then 


V{H = H^ 

The operator A is defined then as follows: 

xi eV(A)nHi 

x 2 eV(y/A) 


V(A) = 


y = 


Xi 

x 2 


A = 


AY = 


0 

—M~ l A 0 

ViX 2 


-M~ l Axi 

thus defined A is closed and has a dense domain. 
Let A* denote the adjoint. Then 

0 -Vi 


A * = 


M~ l A 


0 


= -A. 


To prove this, let 


Y = 


Vi 

V2 


Z 1 

z 2 


; Y,ZeV(A). 
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(4.8) 



Then 


[AY,Z\ e = ^/Ay<2, VAzi -{Ay^z^ 

= [My 2 , M- x Az x \ - [VXy lt y/Az 2 ] 

= [ Y , A'z] e 
= ~[Y,Az} e . 

In what follows we shall omit the subscript E in inner products and norms where 
elements of He are involved. In particular, for Y in the domain of A: 

\AY,Y) + \Y t AY] = Re [AY, Y] = 0. 
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Hence 


||S(f)V'H 2 = llV'll 2 . 

The domain of A being dense, we have that 

l|s(i)y|| = rii, YeH E . 

In particular, the energy stays constant in time. 

More generally, for any Y (0) in He, 

Y (t) = AY (t) + Bu(t) 


interpreted in the weak sense 

x(y(t),y] = [y(0M*y] + [Bu(t),y] 

at 

for every Y in the domain of -A*, has the unique solution (see [Balakrishnan 1981]). 


Y(t) = S(t)Y( 0) + J* S(t-a)Bu(a) da, t> 0 


for 



da < oc 


for every t > 0. 


Remark 

We can finally relate (2.14) to (2.1). Requiring the initial condition vector Y 
to be in V(A) means in particualr that we can define the necessary partial derivatives 
in (2.1) and thus have a solution in the “ ordinary” or pointwise sense. 

Spectral Properties of A 

The spectral properties of A are readily deduced from the spectral properties 
of A. Thus the equation 

AX — AX = Y 

®1 . y 1/1 

’ ’ V2 
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yields for A' : 


Axi - V l x 2 = j/i, 


Xx 2 4- M 1 Ax i = t/ 2 - 


Or, 


) 


\ 2 Mxi 4- Axi = M (y 2 4- Xyi — XVqx 2 ) 

V\x 2 = Ax x - 2/1 

The eigenvalues of .4. are thus given by 

A 2 Mxi 4- 4xi = XVqx 2 
V\x 2 — Axi 

Zero is an eigenvalue of .4., since zero is an eigenvalue of A. In fact 


) 


0 




0 

<t> 


0 


0 


if A<j> = 0. 


Thus the dimension of the null space of A is 6. 
For nonzero A, we note that 


*k = 

is an eigenvector with eigenvalue (to/*). For 


4>k 

iuk4>k 


(4.10) 


(4.11) 


[M<f> k ,V o*] = 0, or <p k eV(A)nHi 
Vo4>k = 0; —u)\M 4>k + A4>k ■- 0, 

and hence (4.11) is satisfied. The dimension of the eigenfunction space is equal to 
the dimension of the eigenfunction space of A corresponding to We shall for 
simplicity take this dimension to be equal to one. Hence the eigenvalues of A may be 
enumerated as: 

w 0 = 0 

{±i|wfc|}, u> k / 0, k > 0; u k -n > u k 

A$k = Wfc 7^0 

wjt / 0, since <f> k = <p k 
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where 


$k = 


iWk\4>k 


4>o ,i. t = 1, . . . ,6 = 


where 

*0,< = *0,i 

[M<t> 0,t, <t>Qj] = 6j- 


A*$k = ~iu k $k 
A*$k = iKI^fc 

A'*o,i = 0. 

These eigenfunctions are orthogonal: 

[**,*,•] = o. 

= 0, 

[$ fc , $*] = <£*] • 

To orthonormalize the we need only to take 

[M<t>k,M = 2^2 ^ 

which we shall assume in what follows. 

Let us show that the {$*,$*:} are complete. Suppose for some ^ in He 

[$*,■0] = [$k,'l> } = 0 for every k. 


Writing 
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we have 


0 = 'I'] = ul[M<t>k,ipi] + iu k [M<pk,i> 2 ]. 

0 = [$*, = vl[M<t>k,Tpi} - iuk[M<t>k,il>2] 

0 = [$o,k,*] = 0 

from which it follows that 

\M4> k ,i> 2 \ = 0 for every A; 
and from the completeness of {<£*} we have that 

ip 2 = 0. 


We also have that 

[M<t>k,^>\\ = 0, Wfc/0 

and by definition is in T-t\. Hence it follows that 


i>i = 0. 


As a consequence we have the modal expansion: 


Y 


00 00 r ® 

+ E [Y,*k]*k + Ei y -^]^o 
1 1 1 


(4.13) 


Since 


S(t)$ k = e^jt; S(t)$* = e-^^; S(i)*o,fc = e^So,* 


we have also 


00 00 

S(t)Y = + £[y,* fc e"***** 

1 1 

+ El*’’ ^°.fc] $o,*- 


(4.14) 


Resolvent of A 
For 


A ^ ±i|u>fc|, A ^ 0, 
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we shall now show that 


XI - A 

has a bounded inverse. To calculate this inverse we go back to (4.8), where we let 


VqX2 = Z- 


Then 


X^Mxi + Axi = Myi + XMyi — XM z. 


For A A ±*|wfc|. A ^ 0, we can define 


w = (X 2 M + A) l {My 2 + XMyi). 


Then 


Hence 


A 2 M (xi — w) + A(xi — w) = -A Mz. 
X^Mixi-Viw) + Aiii-Viw) = X 2 MV 0 w - XMz = 0 


if we let 


z = XVqw 


or, 


x 2 = XV\w -yi + XVqw = Xw - y x . 


Hence we obtain 


X\ 


V\w 

x 2 


Xw - t/i 


(XI - A)' 1 


Vi 

3/2 


where 

w = (A 2 M + A)~ X (M y 2 + Xy{) 


(4.15) 


and 

tu € V(A). 


We use the notation: 


R(X,A) = (A I -A)~ l 
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and refer to the left side as the “resolvent of A .” We have also the modal represen- 
tation: 


KKA)Y = + ^ * 0 , 

^ A — lUk i A + lUk i A 

It is indeed the Laplace Transform: 

R(X,A)Y = r e~ Xt S(t)Y dt, ReA > 0, 

Jo 


(4.16) 


(4.17) 


but of course is defined and analytic in A except for poles at the eigenvalues as follows 
from (4.13) and (4.14). Note in particular that 


B*U{\A)Y = — yi) 


and hence it follows that 

B*n(X,A)B = + A)~ l B = A£* (\ 2 M b + T(A)) _1 B u . (4.18) 

5. CONTROLLABILITY AND STABILIZABILITY 

In this section we show that under a controllability condition, rate feedback 
using a collocated sensor can stabilize the system. All modes will decay even though 
the damping coefficient will decrease with mode frequency. Whatever the initial 
conditions, the elastic energy will eventually dissipate to zero. (This is known as 
“strong” stability.) In our model we neglect any inherent damping in the structure. 
However the controller is robust in the sense it will not destabilize any mode — it will 
only increase the damping. We have always stability enhancement, in other words. 

The main results are known — see [Balakrishnan 1981] and the references 
therein. The presentation here tries as far as possible to be self-contained, using the 
specific features of the problem at hand rather than merely quoting general results. 
In reference to our dynamic equation: 

Y(t) = AY(t) + Bu(t) (5.1) 
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we do not include a damping operator; A is not stable. In fact 


iis(t)y|iE = ni* 

The elastic energy is not dissipated in the absence of control. In finite dimensional 
theory, the system is “exponentially” stabilizable, if (A B) is “exactly controllable 
— by that we mean given any Yi, Y 2 we can find a control «(•) such that 

y 2 = S(t)Yi + f S(t - a)Bu(a) da 
Jo 

for some t > 0. This is impossible in our case because (range of) B is finite- 
dimensional. (See [Balakrishnan 1981], for a proof.) The next best thing we can 
do is to require just “controllability.” Thus we say that (A,B) is “controllable” if 

(J (range of S(t)B) 

t> 0 

is dense in He- This is equivalent to saying that 

IJ ( f S(t - a)Bu(a) da; u(-) e 1/2(0, t)^ 
t>o '*'° ' 

is dense in He • In other words we require that the states “reachable” from the zero 
state are dense in He- 

Theorem 5.1 

A ~ B is controllable in He if ail d only if 

BV ^ 0 (5-2) 


for any eigenfunction <fi, defined by: 

A<j> = oi 2 A/ <(>, 4> 0. 


Equivalently 


^ 0 


for any mode 


>1$ = iu;$, 3> ^ 0. 


(5.3) 
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Proof 


It is convenient to use a modal expansion for any V as: 

Y = f^P k Y + f ]P. k Y + PoY 
1 1 

where P* is the Projector onto the eigenfunction space corresponding to the eigenvalue 
zcjfc, and P_* corresponding to and Po onto the null space of A. Then 

oo oo 

S(t)Y = Y, ei “ ktp * Y + + P ° y - (5.4) 

l l 

If the set 

(J S(t)Bu 

t> o 

where u ranges over all of 7l mc is not dense in He, we can find a nonzero element Y 
in He such that 

[S(t)2?u, Y) = 0, t > 0, u € vr*. 

Hence using (5.4) we must have 

f y^ kt B[u, B*P fc y] + £«“***(«. B*P~ k Y) + [u, B*P 0 Y] = 0, t > 0. 

i i 

But the left side is an almost periodic function in t, and can vanish identically if and 
only if for every k: 

0 = [u, B*P k y] = [u,B'P. k Y] = [u,B*P 0 Y\. 

Since Y is not zero, there must be at least one k such that 

PkY* 0 

or 

PkY = 4 > fc) 

where 

A$ k = iuj k $ k , $>* ^ 0 
[u, B'* k ) = 0. 

But u being arbitrary, we must have 

B*$ k = 0 
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which is a contradiction. Since 4> fc must be of the form: 

<Pk 

$k = 

iuk4>k 

it follows that 

B*4> k = 0 

which is again impossible by assumption. 

Next, suppose A ~ B is controllable. Suppose 


B*4> = 0 


for some mode <fr 

A4> = u) 2 M<p, 4> # 0. 

Then for 

<f> = 

[S(t)Eu,$>] = e^lBu^] = e wt (io;)[u,BV] = 0. 

Hence 3> is orthogonal to 

S{t)Bu 

for every t and u € H mc . Hence (-4 ~ B) is not controllable — contradicting the 
hypothesis. 

We shall show that (A - B) is controllable for our system in Section 4 where 
we study the eigenvalue problem. 



Corollary 

Suppose {A - B) is controllable. Then the number of controls (the dimension 
of the control space) must be at least 6. More generally, the dimension of the control 
space must be at least equal to the largest eigenfunction-space dimension. 
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Proof 


Let <j> 0ti , i — 1, . . . , 6, be a basis for the null space of A. Suppose £?Vo,i are 
linearly dependent. Then 

6 

= 0, not all a* = 0 

l 

or 



But 

6 

&k<t>o,k 

l 

is a nonzero eigenfunction function of A corresponding to the eigenvalue zero and 
controllability is thus violated. 

Next we shall prove the fundamental relationship of controllability to stabiliz- 

abilty. 


Theorem 5.2 

Suppose (.4., B) is controllable. Then the feedback control 

u(t) = a > 0 


is such that 


for every Y (0). 


||y(t)|| — ♦ 0 as t — >• oo 


(5.5) 


Proof 

The closed-loop system dynamic is now 

Y(t) = ( A-aBB*)Y{t ), 
Y (0) given. 


Let 


Y(t) 


Xi(t) 

* 2(0 
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(5.6) 



Then we have 


Xi{t) = V\X<2 (t) 

x 2 (t) — — M~ l Axi(t) — qM _1 BB*xo( 0- 
To relate this to (2.17), let the initial conditions for the latter be given as 


x(0) 

±( 0 ) 


For (3.6) let 


where 


Y( 0) 


xi(0) 

x 2 (0) 


xi(0) = Vix(0) 
i 2 (0) = x(0). 


Define 


Then 


x (t) = xx(t) + Vqx(0) + 



and (5.7), (5.8) yield 


x(t ) = Xi(t) +VoX 2 (t) 
x{t) = x 2 (t) 


(5.7) 

(5.8) 


x 2 (t) = V\X 2 {t) + VoX 2 (t) 

= *(0 

= — M -1 Ax(t) — aM~ l BB*x(t) 


or, 


M x(t) + Ax(t) + aBB'x(t) = 0 


so that the control u(t ) is now 


u(t) = a B*x{t) = -aB*Y(t). 


(5.9) 


(5.10) 
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In other words we have rate feedback using a collocated sensor — a feed- 
back principle for stabilization that is age-old, but still requires proof in our infinite- 
dimensional context. 

First we note that for Y in V(A): 

Re [(>1 — aBB*)Y, V"] = -a\\B'Y\\ 2 

= R e[{A* - aBB*)Y, Y}. 


Hence, see [Balakrishnan 1981], 

(A - aBB *) 

generates a dissipative strongly continuous semigroup. Denoting the latter by 

S a (t), t > 0 


we note that 

l|S«(t)ll < 1. 


Eigenvalues 

Let us consider next the eigenvalues: 

{A-aBB*)Y = A Y 


or, then 
Here 
for if 


(ReA)[y,y] = — a||B*y|| 2 . 

||B*y|| cannot be zero 


B'Y = 0 


(. A - aBB*)Y =AY = XY 

and we violate the controllability condition. Hence it follows that 

-a||£*y || 2 


Re A = 


[Y,Y] 


< 0 . 


(5.11) 
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Hence every eigenvalue has a strictly negative real part. Hence rewriting (5.11) as 

(A I-A + aBB')Y = 0 
and multiplying on the left by 7Z{ A, .4), we have 

y + an(\,A)BB*Y = 0. 

Hence 

B'Y + aB*n{\,A)BB*Y = 0 
or, 

{I + aB*ll(\,A)B)B*Y = 0 

where inside the parentheses on the left side is an m c x m c matrix. Hence the eigen- 
values are the roots of 

D(X,a) = 0 (5-12) 

where 

D{ A, a) = det(7 + aB*H(X, A)B). (5.13) 

But using (4.18), we have 

D{ A, ct) = det [/ + aAB*(A 2 Afj, -I- T(A)) _1 B lt ] . (5.14) 

Let {Afc} denote the eigenvalues, where we know that X k must have the form. 

A k = -| a fcl + iv * > °k,Vk real. 

The corresponding eigenfunctions, denote them Y k , nre then given by 

Y k = H{X k ,A)Bu(X k ) (5.15) 

where 

(I + aM(X k ))u(X k ) = 0, ||u(Ajt)|| = 1 

where 

A* (A) = B'K(X,A)B. 


51 



We note that the dimension of the eigenfunction space is the dimension of the eigen- 
vector space corresponding to the eigenvalue zero of the m c x m c matrix: 

I +aM{\k)) 

and is thus less than m c . It is equal to one, if we assume that the matrix has distinct 
eigenvalues, which we shall, for simplicity, in what follows. 

At this point we leave open whether the sequence {A*} is finite or not. (We 
shall eventually see that it is not.) 

The eigenfunctions { VT} it must be noted are not orthogonal. We shall nor- 
malize them so that 

linii = i. 

Lemma 5.1 

jt,(-0k) < aTrBB* = aT VB*B U . (5.16) 


Proof 

We follow essentially [Gohberg-Krein 1969, p. 101]. Let us orthogonalize { Y)c} 
following the Gram-Schmidt procedure. Then 

k-l 

Zk = Yk — ^ QkjYj) [ %kZj] = 0, k ^ j- 
j = i 

Recall that the {Y*}, k = 1, . . . , n, cannot be linearly dependent for any n, since the 
{A*} are distinct. Hence 

ll^ll#0. 


Now 


Hence 


fc-i 


(. A-aBB*)Z k = A k Y k - ^a fc A^. 

j=i 


[( A-aBB')Z k , Z k ] = A k [Y k ,Z k ] = X k [Z k} Z k \. 


Hence 

D x -a[ BB*Z k , Z k ) 

R k [ z k ,z k I 
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Now 1 SB* is nuclear, B being finite-dimensional. Since 

{ z ‘i ) 

is now an orthonormal sequence, we have: 

“ \BB-Z k ,Z k ) < 

V I 2 *' z ‘] - 

We note that 

TrBB* = Tr B*B = Tr B* U B U . 
Hence (3.9) follows. We have an obvious corollary: 

Corollary 

If the sequence {A^} is not finite, then 

|CT/t| — ► 0 as fc — ► oo 
oo as fc — > oo. 


Proof 

The first part of the statement follows from (3.16) and the second part from 
the fact 

det(7 + a.M(A)) 

is an analytic function for Re A negative, so that no subsequence of the sequence {X k } 
can have a finite limit point. Hence 

|Afc| — ► oo as fc — oo. 

This will also follow from the fact that the resolvent 

U(X, A - aBB*), A ^ X k 

is compact, as we shall see presently. Nevertheless, we still have to prove that the 
sequence { A fc } is not finite, which we shall in Section 7. 
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Now 


S a (t)Y k = e~ |fffc|t e i# ** t yjfc, t > 0 

and hence each mode is damped with damping coefficient <7*. However, the number 
of modes is not finite, and the damping coefficient eventually goes to zero, so that we 
cannot guarantee a finite gain margin. In particular the fact that each mode decays 
is not enough to prove (5.6). 

For this purpose we can invoke a general result due to [Benchimol 1978). We 
are assuming that A — B is controllable and we have seen that the resolvent of A is 
compact, and that 

A + A* = 0. 

Hence by a theorem of Benchimol — [Benchimol 1978] — it follows that the semigroup 
S Q (-) is “strongly stable”: 

||S Q (t)Y|| -» 0 as t -» oo. 

Since 

Y(t) = S a (t)Y(0), 

the theorem is proved. We shall give an independent and self-contained proof in 
Section 7. 

Corollary 

The solution of (5.9) is such that the total energy — elastic plus kinetic: 

W'/A x(t)|| 2 + [ Mx(t ), i(t)] 

is monotonic nonincreasing as t increases and decays to zero as t increases without 
bound. Moreover 


V 0 (*(t» - Po(*(0)-x(0)), as t — ► oc . (5.17) 
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Proof 


We have only to note that 


Y(t) 


V 1 x(t) 

x(t) 


and 

||y(t)|| 2 = \\'/Ax(t)\\ 2 + [Mx(t), x(t)\. 


Thus the elastic energy decays to zero. The rigid-body component is given by 

Voxit) = f Vo x{a) da + V o x(0) 

Jo 

where we so far only know that 


||Po±(t)|| ^ 0 as t — ► oo. 


If the initial state is such that 


is an eigenfunction, then 


vmo) 

i(0) 


Y(t) = e Xkt Y k 


so that 

V 0 x(t ) = \ k e Xkt V o x(0) 


and hence 

Vo(x{t)) = Vo(x(0)) + e Afct P 0 (i(0)) - V 0 (x{0)) 


-> V o (x{0)) - 'Po^CO)) as t — oo. 

We shall show in Section 7 that this holds generally, using the modal expansion. 


Resolvent 

Let us see how the resolvent of (A — aBB*) can be expressed in terms of the 
resolvent of A. We have 

ft(A, A - aBB*)X = Y 
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or 

(XI - A + aBB*)Y = X. 

For A 7^ Afc and A ^ iu>k, 0, we can multiply on the left by Tl(X,A) and obtain 

Y + aTl(X,A)BB'Y = H(X,A)X. 

Hence 

B'Y + aB'R(X,A)BB*Y = B'n(X,A)X 
(I +aM(X))B*Y = B t TZ(X,A)X 
and the matrix on the left side being nonsingular, 

B'Y = (I +aM(\))- l B*R(\,A)X. (5.18) 

Hence it follows that 

Tl( A, A-aBB ') = n(X,A) - all(\,A)B(I + aM(\))~ l B'n(\,A). (5.19) 

It follows in particular that the resolvent of (A — ccBB*) has all the properties of 
TZ(X,A) such as being compact, Hilbert-Schmidt, etc., being a perturbation of the 
latter by a finite-dimensional operator. 


6. ASYMPTOTIC MODES 

In this section we examine in more detail the modes, both open-loop system 
— the undamped structure, and closed-loop — with rate feedback, in particular 
obtaining asymptotic estimates. 

As we have seen in Section 2, the open-loop mode frequencies are the roots of 

det[-u ) 2 Mi + TM] = 0 (6.1) 

where T(X) is self-adjoint, and T(iuj) being a functioi of u 2 . The mode shape is 
determined by 

£(iu k )L(iuj k )- 1 b(iu Jk ) (6.2) 
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where 


(-u) 2 k M b + T(iu) k ))b{iLj k ) = 0. 


Our first step is to show how the dimension of the matrix that determines the 
eigenvalues can be reduced. In (6.1) this dimension is 6m x 6m. We shall show that 
it can be reduced to 6 x 6, regardless of how large m is. 

We shall consider actually the closed loop system, with the corresponding 

eigenvalue problem: 

(A-aBB*)Y = A Y 

proceeding in a slightly different way than before. Let 


Then we have 


Hence 


Aj/i = "P ij/2i A 0 

At/2 "l" A/ -1 A\)2 + ctM 1 BB j/2 = 0 


A 2/2 + t M ^ Ay 2 + 1 B B 1/2 — 0 


or, we need to find 1/2 satisfying 


and then 


A 2 Mj /2 + Ay 2 + ctBB* 2/2 — 0 


V\ = 


We shall call 2/2 the mode shape even though it is not purged of rigid-body modes, 
as (6.5) is. The advantage in going to (6.4) is that we get the undamped mode 
frequencies by setting a = 0, which we cannot do with (5.13). 

To proceed with (6.4), let 


2/2 = 
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We have 


- A 2 f"(s ) + AJ'(s) + Aof(s) + A 2 M 0 / =0, s* < s < s i+1 (6.6) 


X*M h b + a\B u B' u b + A b f = 0. 


(6.7) 


As we have seen in Section 2: 


B u B*b = Db 


where D is the “diagonal” in the sense that 


Db = 


D 1*1 
D m b m 


, Di, 6x6, self-adjoint, nonnegative definite, diagonal. 


f(s) 

ns) 


= e A(\)(s-Si) 


f(Si) 
/'(*« +) 


, Si < S < S i+ i 


where 


0 / 

A>i j (Aq + A 2 Mq) A2 1 Ai 


A(X) = 

The boundary conditions (2.7) relate /'(si+) to / (sj) and /'(s*— ): 


f\ 0+) = AJ 1 (-Li + qAZ?! + A 2 M„,o) /(0) 
/'(*+) = /'(*-) + A % 1 (aXDi + f(si) 


( 6 . 8 ) 


2 < i < m — 1 


and for i = m: 

f\L) = -AZ 1 (la + Aa£ m + A 2 M 6 , L ) /(L). (6.9) 

Hence it follows that we can calculate f(L) and f(L) in terms of /(0) and then invoke 
(6.9) to obtain 

h(\,a)f( 0) = 0 
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where 


h( A; a) = A 2 1 (L\ + ocXD m + A M^i) I 


3 *4.(A)(L— 5 m _i) 


j 4 2 1 (aA£> tn _i + A 2 Mt, m -i) I 


o^(A)(*i+l-*») 


0 


... g- A ( A) *2 

where / is the 6x6 Identity matrix. 
Let us use the notation 


j4 2 1 ( (xXDi + X^Mb,i) I 

I 

j4 2 ^( — L\ ■+■ ot\D\ ■+■ X?Mb,o) 


so that 


Let 


^m(A) O;) 

Then we can write h(X, a) as: 


A i = s i+ i~Si, t = 1, . . . , m — 1 

Ai = s 2 — Sl = S2 
A m _i = L — s m - 1. 

ti(A, a) = Al 1 (A 2 M M + aXDij , 1 < t < m 

/ 0 

U(X,a) I 

Ti(A;a) = | >1 2 1 L 1 + t m (A, a) I | 

I 

— j4 2 1 I»i + tl(A, Of) 


Ti(A,o) = 




fc(A,a) = T m (A;a)g(A > a)T 1 (A;a) 


where 


Q(X,oc) = ^ WA — r „_i( A;Q ) - e« x ^T 2 (\,a) ■ e 


AW Ai 


(6.10) 


( 6 . 11 ) 

( 6 . 12 ) 


( 6 . 13 ) 
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We note that the <*( A; a) as well as Ti{ A; a) are polynomials in A while 


e AW&> 


involve transcendental functions of A. 

The coefficients of the terms of the highest degree in A as well as a are contained 
in the term 

A 2 (aXD m + A A/&,l)Pi2(A! A m - 1 ) • A 2 + A Afj jm _ x ) • ••• 

• P l2 {\-A 2 )A2 1 (<x*D 2 + \ 2 M bi2 ) ■ P 12 (X;A 1 )A2 1 (aXD 1 + X 2 M bfi ). (6.14) 
In particular, the term containing highest powers of A that occurs is: 

A 2m A2 l M btL P 12 ( A; A m _i) A^ v M b , m -x • • • Pi 2 (A; A x ) A^M^. (6.15) 

The term containing the highest powers of a that occurs is 

(x m A2 1 B m P 12 (X; A^A^D^ ■ • • P 12 (A; A^^i) . (6.16) 

From (4.13) we have that 

||A(A;a)|| < ||r 1 (A;a)|| ... ||T m (A;o)|| • ||<^ A|A '|| || e ^' A "-‘|| (6.17) 


where || - || denote matrix norm. 

The corresponding mode shape function / (•) (corresponding to y 2 ) is given by: 

f(s) = | I 0\e AW{s ~ Si) Ti(X-a) • e^^-'T^X; a) 

I 

■■■ e AWAi /( 0), 

A 2 1 (—Li +- A 2 M{ () o + A aD) 

Si < s < Si+ 1- ( 6 . 18 ) 


Let 


d(A;a) = det/i(A;a) 


Then d( A; a) is an entire function of A and the eigenvalues {A*} of (.4. — aBB*) are 
the nonzero roots of 

d{ A; a) = 0. 
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Behavior at A = 0 

Let us first consider 


A = 0. 


We have: 


h( 0,a) = h(0, 0) 


A% 1 Li 


I e AW 


I 

A2 l L\ 


0. 


Moreover we have the power series expansion about zero: 

^2 

h( A; a) = Xh'(0, a) + — h"{0,a) 4- terms of higher order in A 
and correspondingly 

d( A; a) = A 6 a 6 d 6 (a) 4- terms of higher order in A (6.19) 

where 

<* 6 ( 0 ) ± 0 . 

Thus d( A; a) for nonzero a has a zero of order 6 at A = 0, while d( A, 0) has a zero of 
order 12 at A = 0. 


Relation of d( A; a.) to D( A; a) 

Let us examine next the relation of d( A; a) to D( A, a), the latter defined in 
(5.13). Now 


Ad (A) = B*n(X,A)B 

°° 

! X- iu)k 


= E 


^ B' P_ t B &PpB 

1 A 4- iw/t A 


with Pjt as in Section 3 (see also (4.16)) and hence as A goes to zero 

D(A, a) ~ det ^7 


, aB*P 0 B\ 
det ( / + J 


(!)'*• 


- 4- B*P 0 B . 
a 


(6.20) 
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Now the range space of B* PqB is of dimension 6 and hence the null space is of 
dimension 

?/t c 6 


(which is nonnegative by virtue of the controllability assumption!). Hence 

JTlc 6 


det 


- + B*PqB 

a 


X 

Q 


(nonzero constant). 


Hence it follows that 


D( X;a) (!) • (nonzero constant), as A — » 0. 


In a similar way we see from (6.20) that the nonzero poles of D( A, a) are the zeros of 
d( A, 0) to the same order. The dimension of the eigenfunction space of A for A ^ 0 
may be taken to be unity, since we see from (6.18) that the dimension is equal to the 
dimension of the eigenvector space of h( A, 0) corresponding to the zero eigenvalue, 
and if h( A, 0) has distinct eigenvalues, the dimension is equal to one. We can make a 
similar statement for the eigenfunction spaces corresponding to the eigenvalues with 
nonzero imaginary parts of 

A - aBB*. 


Hence it follows that 

d(A; 0) D(X; a) 

is an entire function with zeros coinciding with that of d(A;a). Hence we know that 
we must have: 

d{ A; a) = e qW d(X]0) D(X;a) (6.21) 

where q( A) is an entire function. 


Order of d(A; a) 

Next we shall show that the order of the entire function d(A; a) is less than or 
equal to one. Let 

m(r;a) = max |d(A, a)|. 

|A|=r 
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Now, the determinant of a matrix being the product of the eigenvalues, we have 

6 

|d(A; a)| < (spectral radius of h(A; q) ) 


< || h(X;a) 


From (6.17) we see that 


log||MA;«)ll < X;iog||r i (A ; a)|| + £log Ilex'll- 

j= 1 * =1 

Since the Ti( A; a) are polynomials in A, we need only to consider the order of 

For this purpose, we proceed to evaluate the eigenvalues of .A(A). 

An eigenvector of .A (A), corresponding to the eigenvalue 7(A), must be of the 

tonn vW 

7(A)y(A) 


where 


Hence 


7(A) 2 y( A) = A^(\ 2 M 0 + Ao)y(X) + 7 (A)A 2 1 A 1 y(X). (6.22) 

7 (A) 2 A2i/(A) = A 2 Moy(A) + Aoy(A) + 7 (A)Aiy(A). (6.23) 

02(A) = [A 2 y(A), y(A)]; o x (A) = [Aiy(A), y(A)] 
m( A) = A 2 mo(A) + [Aoy(A), y(A)] 
mo(A) = [Moy(A), y(A)]. 


Then solving the quadratic equation 


7 (A) 2 a 2 (A) = m(A) + 7( A ) a i( A ) 


we have 


7(A) = 2 Q 1 ( Aj ( +fll ^ ± \/ fl i( A ) 2 + 4m ( A ) a 2( A )) ' 
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From (6.22) we have 


'tW\ 


Since 


and 


—J 2/(A) = ^A 2 1 Mq + r-^Ao + y(A). 


02 (A) > (smallest eigenvalue of i42)||y(A)|| i 


|ai(A)| is bounded, 


it follows from (6.24) that 


7(A) 


is bounded. 


Hence 


( 


i4 2 1 A/o + -^2 -^o + ~ ~ ~ ■^l'j — A% ^ Mo 


-(A) 


as A — *• oo. 


Hence “normalizing” the eigenvectors so that 

l|y(A)l| = i, 

every sequence {y(A n )} has a subsequence which converges to one of the eigenvectors 
of 

A^M o 

as |A n | — » oo. We assume (for simplicity) that the eigenvalues of this matrix are 
distinct, which are of course strictly positive. Let fi,, i = 1, . . . , 6 denote these 
eigenvalues and e* the corresponding eigenvector of unit norm. Then for |A| large 
enough we can arrange so that the eigenvalues of .A(A) are 

7, + (A),7r(A), * = 1, .... 6 


'7?<A)' 


- M 


0 as |A| — *• oo 
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with corresponding eigenvectors: 


y»(A) 


y»(A) 

7, + (A)y<(A) 

1 

7r(A)yi(A) 


||y,(A)|| = 1, * = 1, ... 1 6 , j/i(A) -» e*. 


In particular 


MiI/i(A), j/i(A)] — Miei, ei). 

But the ej being real- valued and 

Ai = —L\ -I- L{, 

we have that 


[4i 2/»(A), »i(A)] — 0. 


And as a result, in (6.24) 


and hence 


7 ± (A) — — \/ Q i (A) 2 + 4m(A)fl2(A) — » 0 as |A| — + oo 

7 ±(A) - ±Av7 ~i -» 0 as |A| — » oo. 


Hence we have 


Hence 


e 7fWA 


e ±\y/mA 


1 as | A | —* oo. 


(6.25) 


spectral radius of e"^ A ^ A = max |e ±A '/' I * A , 
for |A| sufficiently large. Again, since y<(A) — *• e< we can find a constant Af such that 

|e^ A || < Mmax|e ±A ^ A |, 0 < A, (6.26) 

for all | A) sufficiently large. Hence it follows that 


max||e^ A)A ll < e rAuMX ^ 

|A|=r 


and hence 


m— 1 


max log||e^ (A)A ’|| ~ r (max L. 


Ui 
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Hence we have that d( A; a) is of order less than or equal to one — or is of “exponential 
type” in the terminology of [Levin 1980]. 

Since 

|Z)(A,q)| — > 1 as | A ] — * oc, 

we note that 

d(A, 0) D(X, a) 

is also of exponential type and hence [Levin 1980, p. 24] we can sharpen (6.21) to: 

d( A, a) = e plW d(A,0)D(\,a) (6.27) 

where Pi(A) is a polynomial of degree one at most. 

Deadbeat Modes 

An eigenvalue which is real is often referred to as a “deadbeat” mode. They 
occur in closed loop only if there are rigid-body modes (zero eigenvalues in the open 
loop). In fact we have: 

Theorem 6.1 

For each a > 0, the number of deadbeat modes is equal to the number of 
(linearly independent) rigid-body modes ( = dimension of the null space of A). 


Proof 

Let life denote the M-orthogonal projection operator projecting H into the 
eigenfunction space corresponding to the eigenvalue of A. Then we have 


x , . <*B*n 0 B . 

H{ A, a) = / + + 


A 2 + 


uji 


We are only ineterested in A < 0. Hence we can write 


(6.28) 


] A] H(X,a) = | A | - aB'UoB - 5 B*U k B (6.29) 

i '' + u k 


where we note 


E 


A 2 + 


B*n k B 


0 as I A 


0 . 
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Hence 


\\\H(\,a) — —aB'UoB as |A| 0. 


Let 


7 *(A,a), k = l,...,m c 
denote the eigenvalues of H( A, a). Then 


|A| 7 fc(A, a) -* (-a) • eigenvalues of as |A| — » 0. 

Now B*n 0 B has exactly 6 nonzero eigenvalues. For, since B is one-to-one, the range 
space of B*IIoB is the same 

{B*4 >} , 4>e null space of A 


and the latter, as we have seen, has dimension 6, by controllability. Hence 


|A| 7 jt(A,a) — < ► -|7fc|> fc=l, . . ,6, as |A| — ♦ 0. 

Hence 

7 ^ (A, a) — -oo, k — 1, . . . , 6, as |A| -» 0. 
On the other hand 

H{ A, a) — ► / as |A| — + oo 

and hence 

7 *.(A, a) > 0 for |A| > A 0 . 

Hence it follows that 


7 it(A,a) = 0 for some A, — | A 0 | < A < 0. 


Hence it follows that 


D( A fc ,a) = 0, A fc < 0, k = 1 ,..., 6 . 
Or, we have exactly 6 deadbeat modes. 
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Remark 


For a = 0, the eigenvalues are the zeros of d(X, 0). We want to consider now 
the limiting case a = oo. For this purpose we consider now 


a +m ^) 


and note that the zeros of 


det 


(^-+.M(A)) (6.30) 

are the same as those of D( A, a). However the form (6.30) allows us to consider the 
case for large a, or a = infinity. The matrix 


b M. (A) — * M{ A) as a — ♦ oo 

Q 


(6.31) 


for each A ^ to;*, and hence (6.30) 


m(A) = detAl(A). 


Hence we define the eigenvalues corresponding to a = -boo as the roots of 


m(A) = 0. 


(6.32) 


We shall show now that the roots are pure-imaginary. For suppose 


,M(A)u = 0. 


Then 


[R(\)Bu, Bu) = 0. 


Let 

7l(X )Bu = y. 


Then 


[R.(X)Bu, Bu] = [y, A I -AY] = X [Y,Y] - [ Y,AY ]. 


Hence 

Re[.4(A)Hu, Bu] = Re(Aj[y,y] 
(B*Tl(X)B is a “positive-real” matrix) 
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and hence 

Re (a) = 0 

or A is pure-imaginary. 

Let us examine the eigenvalue further. We have: 

AB* (x 2 M b + T(A)) B u u = 0. 

We assume that no zero of m( A) is a zero of L(A). In that case we can write 

\B* U L{\) (\ 2 M b L(\) + K{\j) 1 B u u = 0. 

Since 

T(X)* = T( A) 

(l(A) (a 2 M»L(A) + KWY'J = (x'MtLW + KWY'LW 

= L(X) (x ! Mil( X) + K(X ))" ‘ . 


Hence 

AB* (a 2 Mj,L(A)* + K (A)*) ' L(A)*B„u = 0. 

In the special case where 

m c = 6m, B u = Identity 

we have that the eigenvalues corresponding to a = oo are the “clamped” modes. This 
is not true in general, as in fact the example in Section 8 shows. 

The mode “shape” associated with these eigenvalues is given by 

£( A)a(A) 

<t> = 

L{x)a(x) 

where 

d(X) = ^A 2 Mj, + K(X)^j B u u. 


Then letting 


/ = £(A)a(A) 
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we see that 

A 2 M 0 /+p = 0, (6.33) 

where g is defined by (2.6), and 

BlL(X)a(X) = 0, (6.34) 

or, these are inodes in which the control nodes are clamped. We shall characterize 
these modes more precisely below. 


Root Locus 

From Theorem 6.1, it follows that oscillatory modes of the undamped structure 
remain oscillatory for all values of a, however large. The behavior of the set of 
eigenvalues is such that as a increases from zero they migrate from the imaginary 
axis to the left half-plane and then back to the imaginary axis. It is possible to define 
the eigenvalues each as a function of a and show that the real part decreases first 
and then at a critical value of a starts to increase as a increases, going to zero as a 
increases to infinity. We can also show that the critical value increases as the mode 
number increases. The loci describe differential arcs in the complex plane. Thus let 


A fc (0) = iiv k 


where iu>k is a zero of d( A,0). Then we define Afc(a) using the derivatives at a = 0. 
Thus 

0 ) (6 . 35) 


da 


la=0 dA(iw fc ,0) 

where the subscripts denote partial derivatives, and calculate similarly higher order 


derivatives using the identity: 


d(X k (a), a) = 0. 


We can show that it is real and negative, and in particular leading to an approximation 
for ak, the real part, via the Newton formula: 


dq(Wfc, 0) 

d\(iuk, 0 ) 


for small a. 


(6.36) 
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Owing to space limitations we must stop here; and refer to the example in Section 8 
for more. 

Asymptotic Modes 

The modes are the roots of the equation: 

d{\;a) = 0. 

Our interest is not evaluating the roots — which in a given case will be a problem in 
numerical analysis — but rather in their asymptotic behavior as the mode number 
increases without bound. 

Let { Afc} denote a sequence of modes, where we note that 

| Ajfc| -*oo asi-»oo. 

We shall say that the sequence {A*} is asymptotically equivalent if the sequence 

|Ajt — A/t| 

is bounded. In our case we shall show that it actually goes to zero. We call {A fc } 
“asymptotic” modes. Note that the “percent error” 

— — —>0 as fc — *• oo. 

A fc 

Since all mode determination is approximate only, this is clearly the best we can do. 
The zeros of d{ A, a) for each a > 0 are confined to the strip 

- |ct| < Re A < 0, a = sup |<7j| < oo (6.37) 

and from (6.26) we have the important result that in this strip 

pWA.1, i = 

is bounded. Now we can express h( A; a) as: 

2m 

h( A, a) = ^A fc /ifc(A;a) (6.38) 

o 
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where the coefficient matrices hk(-, •) are bounded in the strip (6.37). The coefficient 
of A 2m is given by (from (6.14)): 

WA;o) = i*2 1 M W ;f>i2(A > A m _ 1 ) P 12 (\, (6.39) 


and does not depend on a. The zeros of d(X, a) are those of 



and hence are “asymptotically” those of 


deth 2m (A;a) 


as |A| — » oo. Since Ai and Mb,i are nonsingular we have that 

m-1 

det/i2m(A;a) = detP^A; A*). 
i 

Next we shall show that: 


(6.40) 


Lemma 6.1 


where 


dWA) ■ S i^ur) 

ht(A) - ky/jik \ — 0, |A| — * oo. 


(6.41) 


Proof 

Using 


Vi k(A) 


Vk( A) 


0 

7fc(A)»*(A) 


-7fc(A)j/fc(A) 


27fc(A)yjk(A) 


we have 

Pi 2 (A;A)y fc (A) = | I 0|e^ A | ° 

mA) 


sinh 7 fc(A)A 

= -^jA- WW 
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Hence 


where by (6.25), 

as required. 

Next, let 

Let 


_ JL sinh7fc(A)A 
det Pi2 ( A; A) = 


|7fc(*) - A \Z7Tfc I -» 0, as |A| — * oo 


^oc(A) = 


0 / 
X 2 A2 1 Mo 0 


A s = y A 2 1 Mq 


where the eigenvalues of j4 s are 


Then 


0 A»WA _ 


y/jik > k = 1, ...,6. 

cosh(Ai4 a A) (Ai4 a A) _1 sinh(AA s A) 

(Ai4 s A) sinh(AA a A) cosh(A-A 4 A) 


Now 


l|Pi2(A;A)|| 

is bounded in the strip (6.37) and hence it follows that for every k 
||(P 12 (A;A) - (A/l Jt A) _1 sinh(A>l 4 A))e/k|| -> 0 


and hence 


||p 12 (A;A) - (Aj4 s A) 1 sinh(Aj4,A)| -» 0 


as |A| —* oo in the strip (6.37). 
The zeros of 

are given by 


sinh 7*(A)A 

7fc(*)A 

7 fc (A n )A = im r. 


(6.42) 
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Let 


X' n ^/^^k A = inn 


being one sequence of zeros of 


det ((A,4 S A) -1 sinh(AA s A)^ 
corresponding to the eigenvalue 

X'nVJik and eigenvector e* . 


Then 


|(7fc(^n) - A n> //Ifc ) A| = |(A n — X' n )y/]Ik A| -» 0, as n — + oo 


and hence 


| A n — A^| — ► 0, as n — ► oc. 


We say in this case { A n } is asymptotically equivalent to { Aj,} — or that the asymptotic 
zeros of det Pi 2 (A; A) are given by 


* ±inn 

*n,k = a > K = 1, . . . , b. 

VMfc A 

Thus the asymptotic zeros of (6.40) are given by 

±inn 


(6.43) 




, fc = 1, . . . ,6, j = 1, . . . ,m — 1. 


VJ*k Aj 

Let us turn now to the asymptotic zeros of rf( A; a). Let 

<521 (A; A) = AP 21 (A;A); Q 12 (A;A) = AP 12 (A;A) 


(6.44) 


and 


h( A;o) = A Tn (q m (A;a) + 


r(A 


r 1 ) 


(6.45) 


where, using (6.14), 


9m(A;o) — -I- Qi 2 (A;A T n— i) 

<5i2(A; Ai) A 2 1 ^Mb,o + — 


(6.46) 
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and 


lk(*;a)ll 


is bounded in the strip (6.37). Hence we can write 

6 1 

d{ A;a) = detg m (A;a) 4- Ylji; Q ) 


where 


d k ( A; a) are bounded in the strip. (6-47) 

For large |A| we can use the approximation (6.42) and hence 

detq m (A;a) = JJdet A% 1 ^Mb,L + II det j(A a Afc) sinh Aj4 s Afcj 


which for all |A| sufficiently large can be expressed: 

(constant) JJ (jj + J ^i( A ! a ) + d *( A; a )) ' ^ 6 ' 48 ^ 

Given e > 0, we can make |A| large enough so that the second term is less than e. 
Hence taking 


7 fc (A)A fc = im + 9, 


where 

| sinhd| < e 


or, approximately 

M17T + 6 

^ Aky/Tlk 

We see that there is a value of 9 such that (6.48) is zero. Since we know that the real 
part of the eigenvalue must be negative, we have that 


Red < 0. 


We see however that the zeros are again asymptotically the same as that given by 
(6.35). In summary asymptotically the zeros of d( A, a) for any a are given by the 
zeros of 

det ^ JJ (A^s ) -1 sinh Aj4 s A^ . 
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Mode Shape 

The (unpurged!) asymptotic mode shapes are determined by the eigenvector 
V- 

h^n{\\a)y -- 0 

where y can be determined as follows. Fix j, and let 

A * n7r 

Am- ly/JPj 

and assuming for simplicity that the A * are distinct, we see that Pi 2 (A n ; A,) are 
nonsingular for i not equal to m — 1, and hence we can take 

V = (^2 Mb,m— 2-Fi2(A n ; A m _2) Pilin', Ai)A 2 Mb,oj 6 j 

so that 

Pi2(^n> A m _i)ej = 0. 

We can repeat this procedure for A m _ 2 , . . . , Aj. The corresponding mode shape is 
then determined by (6.18) where /(0) is now denoted y. 


7. MODAL EXPANSION 

For the undamped structure (a = 0), we have seen that we have a “modal 
expansion” in terms of the eigenfunctions of A given by (4.16). The eigenfunctions 
are orthogonal and complete. The question arises as to what extent this property 
holds in the closed-loop case — for the eigenfunctions of 

A - aBB *, a > 0. 

We have seen that these functions are not orthogonal and hence we need to examine 
what happens to the modal expansion. 

We have seen (cf. (5.15)) that the eigenfunctions are of the form 

Yk = Tl{\ k) A)B'u{\ k ) (7.1) 
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where {A*.} are the eigenvalues. For simplicity we shall assume that the dimension of 
the eigenfunction space is unity. Let Z k denote the eigenfunction of 

(A-aBB*)* = A'-aBB' 

corresponding to the eigenvalue A* . Then the pertinent properties of A * being similar 
to those of A, we have: 

Z k = TZ (a k, >1*) B*v (A*) , u(-) € E™* 

= -n (-A fc , A) B*v (a*) . ( 7 . 2 ) 

The main feature of these eigenfunctions is that they are “biorthogonal” 

[Y k , Zj] = 0, k A 3 
[Y k> Z k ] A 0. 


Riesz Basis 

Recall now that a sequence {$ k } of elements in a Hilbert space H is called a 
“basis” if every element Y in the space can be expressed as 

OO 

y = J^a k V k 

l 

where 

OO 

]£l a *l IM < oo 

i 

and 

OO 

0 = ^Qjfc^fc 

i 

implies 

a k = 0 for every k. 

A biorthogonal sequence {$*, ip k } is called a Riesz basis if there is a linear bounded 
operator T on H into H with bounded inverse such that {T<p k } is an orthonormal 
basis. This implies in particular that we have the expansion: 

y = £ [y = EiyM** (7.3) 

l l 
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where 


£|[ F ,^]| 2 < oo ; 


Elt^H 2 < oo 


T<p k = T * V*. 


Also: For Y, Z in H 


[y,z] = 


(7.4) 


The matin result in this section is that {Y k , Z k } upon “normalization” so that 


[Y k ,Z k \ = 1 


( “normahzing” u(Ajt), v [\ k j appropriately, which we shall assume from now on) form 
a Riesz basis. This will follow from [Balakrishnan 1996] upon verifying the conditions 

i) 


| A* + A j 


> 6 > 0 , 


for k f j 


ii} 

< oo 

1 

iii) (algebraic) multiplicity of each eigenvalue is equal to unity. 

The condition (i) follows readily from the asymptotic estimate (5.35). Condition (ii) 
is verified in Section 5, (5.16). Condition (iii) is automatic since we are assuming the 
dimension of the eigenfunction space for each a > 0 is unity. 

We can now proceed to exploit the modal expansion (7.3). First we note that 
if A*, is an eigenvalue so is A* if 

Im A fc 0 


and there are exactly six real-valued A*. Let 


(A - aBB*)Y k = A k Y k . 


Then 

(A-aBB')Y k - A k Y k 
(. A-aBB')*Z k = A k Z k . 
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where {A*.} are the eigenvalues. For simplicity we shall assume that the dimension of 
the eigenfunction space is unity. Let Z * denote the eigenfunction of 

(A-aBB'Y = A*-aBB* 

corresponding to the eigenvalue Ajt . Then the pertinent properties of A being similar 
to those of A y we have: 

Z k = n (Ifc, >1*) B*V (I*) , «(■) 6 E m ' 

= -n (-At, a) b*v (a*) . (7.2) 

The main feature of these eigenfunctions is that they are “biorthogonal” 

[Y k ,Zj] = 0, k^j 

[Y k ,Z k ] / 0. 


Riesz Basis 

Recall now that a sequence {'iftk} of elements in a Hilbert space l~t is called a 
“basis” if every element Y in the space can be expressed as 

00 

Y = 

l 

where 

00 

El°*l IhMI < 00 

i 

and 

OO 

0 = 

1 

implies 

dk — 0 for every k. 

A biorthogonal sequence ifci t} I s called a Riesz basis if there is a linear bounded 

operator T on H into H with bounded inverse such that {T(j) k } is an orthonormal 
basis. This implies in particular that we have the expansion: 

°o oo 

Y = ( 7 - 3 ) 

1 1 
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where 


00 00 

J 2\[ Y , rP k )\ 2 < 00 ; £|[ r ,<* fc ]| 2 < 00 

1 1 

T<p k = T*~ 1 ip k . 

Also: For Y, Z in H 

OO 

[Y,Z] = E [Y,*k][*k.Z]' (7-4) 

1 

The main result in this section is that {V*, Z k } upon “normalization” so that 

[Y k ,Z k ] = 1 

( normalizing” u(X k ), v (^X k ^j appropriately, which we shall assume from now on) form 
a Riesz basis. This will follow from [Balakrishnan 1996] upon verifying the conditions 

i) 

|a* + Aj > 6 > 0, for k ^ j 

ii) 

OO 

EM < 00 
1 

iii) (algebraic) multiplicity of each eigenvalue is equal to unity. 

The condition (i) follows readily from the asymptotic estimate (5.35). Condition (ii) 
is verified in Section 5, (5.16). Condition (iii) is automatic since we are assuming the 
dimension of the eigenfunction space for each a > 0 is unity. 

We can now proceed to exploit the modal expansion (7.3). First we note that 
if A*, is an eigenvalue so is A* if 

Im X k ^ 0 

and there are exactly six real-valued A*. Let 

(A - aBB')Y k = A k Y k . 

Then 

C A-aBB')Y k = X k Y k 

{A — aBB*Y~Z k = X k Z k . 
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Hence numbering so that Aj, . . . , Ag are real, and the A/~ are 
|Afc|, we can express the modal expansion as: 

y = £iy,z t ]n + Y.([Y,z t \Y t + [y,z t ] 

1 7 

Correspondingly the solution of the closed-loop system: 

Y(t ) = (>l-aBB*)y(t); r(0) = y 

can be expressed: 

y(t) = + ^([y.Ztie^'yfc + [y, z t ] e 4 *' y*) . (7.6) 

1 ^ 

We can easily deduce strong stability of the semigroup S a (-) from (7.6), or 
equivalently from the fact that the eigenfunctions {Y*} are complete in He, exploiting 
the dissipativity: 

I|5.(«)||<1. 

If Y is real-valued, we note that 

Wm = [y,z»] 


in increasing order in 

- n) . (75) 


and hence 


[Y,Z k ]e Xkt Y k + [Y,Z k ]e Xkt Y k 

can be expressed 

= [[Y, Z k ]Y k + [y,Z fc ] Y k ] cos u k t + e- |<Tfc|t [\Y, Z k ]Y k 


[y, z k ] y*] sino/*t. 

(7.7) 


Since we know that we must have: 


Y k = 


Vk 

^kVk 


; = 


Zk 

AjkZfc 


we can proceed to develop expansions for ii(t), x 2 (t) where 


going back to (5.7), (5.8). We omit the details. Prom (7.6) we can readily deduce 
(5.17). 


8. ILLUSTRATIVE EXAMPLE 

To illustrate the foregoing theory and concepts, we consider now an example 
— simplified in the extreme to reduce notational complexity and wholly non-numeric 
to avoid computer calculation. Thus we consider beam torsion about a single axis 
with a control at one end and a lumped mass at the other, and no interior nodes. 

Retaining the nomenclature of Section 2 as much as possible but using 9(t , s ) 
in place of 0(f, s), the dynamics can be described by: 

m^O — cgg0" = 0; 0 < s < L; 0 < t (8.1) 

yielding in the notation of Section 2: 

Mo = 77144; ^2 — C66- 

The abstract version becomes: 

H = L 2 (0,L) x E 2 

/(•)eL 2 (0,L), beE 2 , b 
The stiffness operator A is then given by: 

6 i 2 ( 0,L); 

Ax = y\ y = 9 

c 

<?(s) = egg g (s), 0 < s < j, 

-c 66 /'(0) 
c 66 /'(L) 



bo 

b L ' 
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Thus defined 


[Ax,x] 


-cee [ L f"(s)f(s)ds + c 66 f(L)f"(L ) - o 66 /(0)/'(0) 

Jo 

c 66 f L l/'Wf ds (8.2) 

Jo 


yielding the potential energy, as required. There is a rigid-body mode: 


Ax = 0, 


where x is of the form 


x — 


(8.3) 


where 


f(s ) = a, 0 < s < L. 

Placing the control at s = 0, the control operator B is given by 

0 

B u u 


Bu = 


; B u u = 


For x in T)(A), 


x = 


we see that 

B*x = /( 0). 

Since there is a control at one end, we see that all modes are controllable and that 
(A rsj B) is controllable. 

Finally, the mass operator M is given by 

mo 0 
0 m L 


Mx — y; y = 


m 44 

M„b 


; M b = 


Correspondingly we have the boundary equations: 

mo6>(£,0) — C(, b 0'(t, 0) + u(£) = 0 

m L 6(t,L) + ceeO'(t,L) = 0 


(8.1a) 
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The space Hi (M-orthogonal to the null space of A) consists of elements of 
the form 

/(•) 


x = 


bo 

b L 


where 


m 44 f(s) ds -t- m 0 bo 4- m L bL = 0. (8.4) 

The domain of \/~A , by [Balakrishnan 1990] is characterized by elements of the form 


i = 


/(•) 
/( 0 ) 
f(L) 


where /(•) is absolutely continuous and /'(•) G £ 2 ( 0 , L) and in particular 

llv/Ixll 2 


is the potential energy given by (8.2). Thus 

X>(\/7)n?< 1 = 


/ 

/( 0 ) 

/(£) 


/ is absolutely continuous with /' 
in 1/2(0, L), and (8.4) holds 


Also 


H e = (2?(v/I)n Hi) x H 


with energy inner product, as in Section 2. 
The feedback control is: 


u(t) = aB*x(t) = ctQ(t) 0), q > 0. 


Closed-Loop Modes 

We proceed directly to characterize the closed-loop modes. In the notation of 


Section 6: 


0 


1 



where 


r> 77^44 

V = 


- 


cosh Xus 

A^sinhAi/s cosh Xu s 


where 


We see that 


A 2 mr 

1 e A W L 

1 

C66 

aA+A 2 mo 



C66 


MA.a) = 1 I e ' ’ 

C66 

= A (ai(A) sinh XuL + 02 (A) cosh XuL) 


i(A) = —V + Aam L + A 2 m 0 mi,) 
"Cfo V 

a + A(m L + m 0 ) 

02 (A) = 

C66 


d( A; a) = h{X;a) 

is an entire function of order one. It has a zero of order one at A = 0, for nonzero a, 
and of order two for a = 0. It is of “completely regular growth” in the terminology 
of [Levin 1980]: 

log |d(re* e , a) | 

lim = vL\ cos0| 

r — ► 00 7 * 


and hence [Levin 1980, p. 169]: 

N(r) 1 f** 


lim 

r — >0 t 


= -f 

2n Jo 


uL\cos9\ d6 > 0 


where N (r) is the number of zeros in the circle of radius r. Hence the number of 
zeros is not finite. For large |A|: 

d( A; a) ~ X z mQmi sinh XuL 

and hence the asymptotic modes are the roots of 

sinh XuL = 0, 
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or 


±m7r 

~Tl~ 


A 


( 8 . 6 ) 


for all a > 0 . Here however we can make a more exact calculation. Thus the 
eigenvalues {A^} are the roots of 


tanh A i/L + 6(A;a) = 0 


( 8 . 7 ) 


where 


and 


b(X-a) = 


We can rewrite ( 8 . 4 ) as 


77144 a + H m L + m o) 
v 77744066 + Aamx, -|- A 2 mom£ 

|6(A;a)| = O as |A| — 00. 

\vL + tanh -1 i>(A;a) = 0 


where 


and since 


*(A;cO = 


02(A) 


oi(A) 

tanhx = tanh(x ± 2tn7r), n integer 


( 8 . 8 ) 


we have 

XuL = ±m7T + 3 log | + !m ,Q | ~ °> ( 8 . 7 a) 

* l-6(A;a) 

using the principal value of log x, real when x is positive. For | A| large, ( 8 . 7 a) becomes 

XvL = itnir + \ ( 8 . 9 ) 

v \ rriLTno / A 

yielding a slightly better approximation than (8.6), for large n. For nonzero a, we 
can see that ( 8 . 4 ) has exactly one real root, approximately: 

A = . ~ Q - • ( 8 . 10 ) 

777 44L + mi + 777 0 
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Clamped Modes 

We can calculate that 


L(\) = 


1 


cosh XvL 

and the clamped modes are the zeros of sinh XvL or, 

±ikn 


0 

sinh A i/L 
Xu 


\ k = 


vL 


k = 1 , 2 ,.. 


( 8 . 11 ) 


Dynamic Stiffness Matrix 

We can calculate that 


K{ A) = C66 


0 


-1 


Xv sinh A vL cosh XvL 


and hence that 


T( A) = 


Al>C 6 6 


cosh XvL 
-1 


-1 

cosh A vL 


sinh XvL 

which is clearly nonnegative definite for A real, and nonsingular except for A 
Also the inverse of the dynamic stiffness matrix: 


= 0. 


(A 2 mj, + T( A)) 


-l 


A/i(A; 0) 
Hence 


A 2 mi,sinh XvL 4- cgeAt'cosh XvL — Xvcee 

— Xvcm X 2 mo sinh XvL + C66-W cosh XvL 


AS; (a 2 M 6 + T(A)) 1 Bu = — ^t-(a 2 m£,sinhAi/L + cgeA^ cosh XvL} 


h{ A; 0) 
= m(A). 

Hence we can verify that for this example 


(8.12) 


In other words in (6.27) 


h{ A; 0) = h(X, 0) D(A; a). 


Pi( A) = 0 
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which we may conjecture holds in general. 


Root Locus 

Beginning first with the limiting eigenvalues as a goes to infinity, given by the 
roots of 

m( A) = 0; 

we have from (8.9) 

tanhAi^L 4- ^° 66 = 0 (8.13) 

attil 

or directly from (8.7) by taking the limit as a goes to infinity in (8.8). The roots are 
of course pure imaginary: 

An = ipn, Pn real, 

0 n vL = ±inn 4 - i6 n , |fi n | < n. 

These are the modes which satisfy 

m L 9(t,L ) 4- cse9'{t, L) = 0; 0(t,O) = O 

m^9(t,s) - c^9"{t,s) = 0, 0 < s < L. 

These are not the clamped modes, although they are, asymptotically. 

Since we are only interested in the nonzero eigenvalues, let 

F(A;a) = (t /2 Cg 6 4- A ami 4- sinh XvL 4- vc^{a 4- A(m£, 4- mo)) cosh XvL 

whose zeros are the nonzero eigenvalues. Let {«<;*} denote the zeros for a = 0. Fix 
k. Now 

F(A(a); a) = 0 (8.14) 

defines an implicit function A* (a), with 

Afc(0) = iu> k 


and we define all derivatives at a = 0 using (8.14). In particular 


dXkja) 

da 


Q=0 


— F a (Afc(q); a) 
F A (A*(a); a) 


a=0 
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where the subscripts denote partial derivatives, and the main point is that it is real, 
and it is negative. In particular this shows that the real part is decreasing. Since we 
know that the real part goes to zero as a goes to infinity, we see that there is a value 
of a at which its derivative must change sign. 


Closed-Loop Mode Shapes 

Following Section 7, the (unpurged) closed-loop mode shape corresponding to 
the eigenvalue A* is given by 


/*(*) = |1 Ole***" 


l 

aAfc-hA^mo 

C66 


(coshAfci's + ( Q + ^ m o) /(0). 

V "C 6 6 / 


Since arbitrary multiplicative constants can be used, we may define the mode shape 
as: 

fk(s) = Ak s\nh(X k^s + $k), 0 <s<L 


where 


tanh 9 k — 


^66 


a + Afcm 0 

For Xk 7 ^ 0, the purged version would be in the notation of Section 7, 


Y k = 


4>k 


4>k 


; Z k = 


Xk4>k 


Xk<t>k 


where (new notation, not to be confused with Section 2): 


4>k 


hi') 
hi 0 ) 
hiD 


his) = fkis ) - ^TTiu fkis) ds + mofkiO ) + m L f k iL ) 

The constant Ak can be determined to normalize the biorthogonal system as in Sec- 
tion 7 and thus obtain a Riesz basis for He- 
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Limiting Case: m L = oo 

We illustrate finally how to handle the case when one end is clamped. We set 
ttil, = +oo. This results in the boundary condition 

9(t,L) = 0 

replacing the condition at L in (8.1a). We may take 

n = L 2 (0,L) x E 1 


-<W"0 


/(■) 

-C66/'( 0) 

, X — 

/(0) 


and 

Ax = y\ y = 

There are no rigid body modes and the eigenvalues are roots of 

(a + Amo) sinh \vL + vc§§ cosh \vL = 0 


(8.15) 


or 


A n j/X = ±in7r - tanh x 6(A n ) 


The eigenfunctions are 


6(A) 


^C66 

Q! + Amo 


Y k 


4>k 

X k 4>k 


4>k 


/*<•) 

A(0) 


fk(s) = A k sinhX k v(L - s). 


The root-locus problem becomes much simpler than before. Again, we omit the 
details. For a = oo, the modes are the zeros of sinh XvL or 


From (8.15) we see that 


±in 7r 
vL ' 


d A n = -1 

da i _ ( ^c 6 6 "\ 2 

\a+X n mo J 
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and is real negative at a = 0 and goes to (—1) as a goes to infinity. For large n we 
have the approximation: 

“ '"n < o 

— ~ -(^ 66 ) 

and 


/ \2 ^ 


da 


(a 4- A„mo) 2 


Re = 0 

da 


implies 


a c rit = mo(|ReA n | + |ImA n |) 

and shows that the critical value of a increases with mode number. 
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NOMENCLATURE 


A 

B 

s 

(•£ 1 j X 2 , 2-3) 

U 

V 

w 

<h 
4 > 2 

4 > 3 

M 

M 0 

M b 

n 

n x 

L 2 ( 0 ,L) e 
R 6 , E 6 
m c 

[. 1 
V(A) 

/0, /e, 
u(t) 

Ait 

Ok 

Ulk 

Tr 
Re 2 
Im 2 
det M 
|z| 

l|F|| 

1 

A* 

SCOLE 


stiffness operator 
control operator 
position along beam 
rectangular coordinates 
axial displacement (xi -component) 
displacement (xj-component) 
displacement (i3-component) 
torsion angle about Xi-axis 
torsion angle about X2-axis 
torsion angle about X3-axis 
mass/inertia operator 
mass/inertia matrix 

composite matrix of mass/inertia at nodes 
Hilbert space 

space M-orthogonal to null space of A 

L2-space of 6 x 1 vector functions over (0, L ) 

Euclidean 6-space 

number of control inputs 

inner product 

domain of operator A 

moments of inertia 

control input 

eigenvalues 

Real part of A* 

angular mode frequencies 

trace 

real part of 2 
imaginary part of 2 
determinant of M 
absolute value of 2 

norm of vector F ; operator norm of matrix F 
conjugate of 2 
adjoint of A 
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